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Abstract 

These lectures present a survey of recent developments in the area of ran- 
dom matrices (finite and infinite) and random permutations. These probabilistic 
problems suggest matrix integrals (or Fredholm determinants), which arise very 
naturally as integrals over the tangent space to symmetric spaces, as integrals 
over groups and finally as integrals over symmetric spaces. An important part 
of these lectures is devoted to showing that these matrix integrals, upon apro- 
priately adding time-parameters, are natural tau-functions for integrable lattices, 
like the Toda, Pfaff and Toeplitz lattices, but also for integrable PDE's, like the 
KdV equation. These matrix integrals or Fredholm determinants also satisfy Vi- 
rasoro constraints, which combined with the integrable equations lead to (partial) 
differential equations for the original probabilities. 
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introduction 

The purpose of these lectures is to give a survey of recent interactions between statistical 
questions and integrable theory. Two types of questions will be tackled here: 

(i) Consider a random ensemble of matrices, with certain symmetry conditions to 
guarantee the reality of the spectrum and subjected to a given statistics. What is the 
probability that all its eigenvalues belong to a given subset E ? What happens, when 
the size of the matrices gets very large ? The probabilities here are functions of the 
boundary points q of E. 

(ii) What is the statistics of the length of the largest increasing sequence in a random 
permutation, assuming each permutation is equally probable ? Here, one considers 
generating functions (over the size of the permutations) for the probability distributions, 
depending on the variable x. 

The main emphasis of these lectures is to show that integrable theory serves as a 
useful tool for finding equations satisfied by these functions of x, and conversely the 
probabilities point the way to new integrable systems. 

These questions are all related to integrals over spaces of matrices. Such spaces can 
be classical Lie groups or algebras, symmetric spaces or their tangent spaces. In infinite- 
dimensional situations, the " cxD-fold" integrals get replaced by Fredholm determinants. 

During the last decade, astonishing discoveries have been made in a variety of di- 
rections. A first striking feature is that these probabilities are all related to Painleve 
equations or interesting generalizations. In this way, new and unusual distributions 
have entered the statistical world. 

Another feature is that each of these problems is related to some integrable hier- 
archy. Indeed, by inserting an infinite set of time variables ti,t2,t3, ... in the integrals 
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or Predholm determinants - e.g., by introducing appropriate exponentials e^i in 

the integral - this probability, as a function of ti,t2,^3, satisfies an intcgrable hier- 
archy. Korteweg-de Vries, KP, Toda lattice equations are only a few examples of such 
integrable equations. 

Typically integrable systems can be viewed as isospectral deformations of differential 
or difference operators L. Perhaps, one of the most starthng discoveries of integrable 
theory is that C, can be expressed in terms of a single "r-function" r(ti, t2, •••) (or vector 
of r-functions), which satisfy an infinite set of non-linear equations, encapsulated in a 
single bilinear identity^ . The ti account for the commuting flows of this integrable 
hierarchy. In this way, many interesting classical functions live under the same hat: 
characters of representations, ©-functions of algebraic geometry, hypergeometric func- 
tions, certain integrals over classical Lie algebras or groups, Predholm determinants, 
arising in statistical mechanics, in scattering and random matrix theory! They are all 
special instances of ''r-functions" . 

The point is that the probabilities or generating functions above, as functions of 
ti,t2, ■■■ (after some minor renormalization) are precisely such r-functions for the cor- 
responding integrable hierarchy and thus automatically satisfy a large set of equations. 

These probabilities are very special r-functions: they happen to be a solution of yet 
another hierarchy of (linear) equations in the variables ti and the boundary points q, 
namely r(t; c) = 0, where the form -roughly speaking- a Virasoro-like algebra: 

The point is that each intcgrable hierarchy has a natural "vertex operator^^ , which auto- 
matically leads to a natural Virasoro algebra. Then, eliminating the partial derivatives 
in t from the two hierarchy of equations, the integrable and the Virasoro hierarchies, 
and finally setting t — 0, lead to PDE's or ODE's satisfied by the probabihties. 

In the table below, we give an overview of the different problems, discussed in this 
lecture, the relevant integrals in the second column and the different hierarchies satisfied 
by the integrals. To fix notation. Tie, St, % refer to the Hermitian, symmetric and 
symplectic ensembles, populated respectively hj i x i Hermitian matrices , symmetric 
matrices and self-dual Hermitian matrices, with quaternionic entries. Ti.i{E), Sj^{E), 
%{E) are the corresponding set of matrices, with all spectral points belonging to E. 
U{£) and 0{i) are the unitary and orthogonal groups respectively. In the table below, 
Vt{z) := Vo{z) + '^tiz\ where Vo{z) stands for the unperturbed problem; in the last 
integral Vt{z) is a more complicated function of ^1,^2, ■■■ and z, to be specified later. 



4 



§1. p.5 



Probability problem 


underlying t-perturbed 
integral, r-function of — > 


corresponding 

intcgrable 
hierarchies 


P(M G HniE)) 




Toda lattice 
KP hierarchy 


P{M G SniE)) 




Pfaff lattice 
Pfaff-KP hierarchy 


P{M e %{E)) 




Pfaff lattice 
Pfaff-KP hierarchy 


P((M^ Mo) e 

HniEi) X nn{E2)) 


g- Tr{Vt {Mi)-Vs (M2)-cMiM2) 


2d- Toda lattice 
KP-hierarchy 


P{M e HooiE)) 


dct {I-Kt{y,z)lE4z)) 
( Frcdholm determinant) 


KdV equation 


longest increasing sequence 
in random permutations 




Toeplitz lattice 
2d- Toda lattice 


longest increasing sequence 
in random involutions 




Toda lattice 
KP-hierarchy 



Acknowledgment: These lectures represent joint work especially with (but also in- 
spired by) Mark Adler, Taka Shiota and Emil Horozov. Thanks also for many informa- 
tive discussions with Jinho Baik, Pavel Bleher, Edward Frenkel, Alberto Griinbaum, 
Alexander Its, espacially Craig Tracy and Harold Widom, and with other participants 
in the semester at MSRI. I wish to thank Pavel Bleher, David Eisenbud and Alexander 
Its for organizing a truly stimulating and enjoyable semester at MSRI. 

1 Matrix integrals, random matrices and permuta- 
tions 
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1.1 Tangent space to symmetric spaces and associated random 
matrix ensembles 

Random matrices provided a model for excitation spectra of heavy nuclei at high ex- 
citations (Wigner [7^], Dyson |27] and Mehta based on the nuclear experimental 



data by Porter and Rosenzweig [Q; they observed that the occurrence of two levels, 
close to each other, is a rare event (level repulsion), showing that the spacing is not 
Poissonian, as one might expect from a naive point of view. 

Random matrix ideas play an increasingly prominent role in mathematics: not only 
have they come up in the spacings of the zeroes of the Riemann zeta function, but 
their relevance has been observed in the chaotic Sinai billiard and, more generally, in 
chaotic geodesic flows. Chaos seems to lead to the "spectral rigidity", typical of the 
spectral distributions of random matrices, whereas the spectrum of an integrable system 
is random (Poisson)! (e.g., see Odlyzko and Sarnak []60| ). 



All these problems have led to three very natural random matrix ensembles: Hermi- 
tian, symmetric and symplectic ensembles. The purpose of this section is to show that 
these three examples appear very naturally as tangent spaces to symmetric spaces. 

A symmetric space G/K is given by a semi-simple Lie group G and a Lie group 
involution a : G ^ G such that 

K = {x e G, a{x) = x}. 

Then the following identification holds: 

G/K = {ga{g)-' with g e G}, 

and the involution a induces a map of the Lie algebra 

o"* : — ^ such that {cr^Y = 1, 

with 

Q = i®P with I ^ = ^ S s^c^ ^*(^) = ^} 
l^p = {aG0 such that a* (a) = —a} 

and 

[fi, t] c [t, p] c p, [p, p] c t 

Then K acts on p by conjugation: kpk~^ C p for all k E K and p is the tangent space 
to G/K at the identity. The action of on p induces a root space decomposition, with 
a being a maximal abelian subalgebra in p: 



p = a + pa, with nia = dimpQ,. 



Then, according to Helgason [35 1, the volume element on p is given by 
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where A+ is the set of positive roots; see ^ This will subsequently be 

worked out for the three so-called A„-symmetric spaces. See also Sarnak's MSRI-lecture 
||59| in these proceedings, who deals with more general symmetric spaces. I like to thank 
Chuu-Lian Terng for very helpful conversations on these matters. 

Examples : 

(i) Hermitian ensemble 

Consider the non-compact symmetric spac^ SL{n,C)/ SU{n) with a{g) = g~^~^. Then 

SL{n,C)/SU{n) = {gf \ g G SL{n,C)} 

= {positive definite matrices with det = 1} 

with 

K = {ge SLin, C) | aig) = g} = {g E SL{n, C) \ g-^ = f} = SU{n). 

Then a^{a) = —aJ and the tangent space to G/K is then given by the space p = Hn 
of Hermitian matrices 

s/(n, C) = 6 © p = su{n) © Hn, i.e., a = oi + 02, ai G su{n), 02 G Hn- 

If M G Tin, then the Ma, ^Mij and '^Mij {I < i < j < n) are free variables, so that 
Haar measure on M G Tin takes on the following form: 

n 

dM:=Y[dMii Y[ [d'^Mijd'^Mij). (1.1.1) 

1 i<«<i<" 

A maximal abelian subalgebra a C p = 7i„ is given by real diagonal matrices z = 
diag(2;i, . . . , Zn)- Each M G p = 7i„ can be written as 

M = e^z e"^, eK = SU{n), 

withg 

A= ^ {ak£{ek£ - eik) + ibk£{ek£ + Cik)) e t = su{n), a« = 0. (1.1.2) 

l<k<e<n 

Notice that — cek and i{eke + e^) G I = su{n) and that 

- eek, z] = {zi - Zk){ek£ + e^) E p = Hn 



^The corresponding compact symmetric space is given by {SU{n) x SU{n))/ SU{n) 
"^Ckt is the nx n matrix with aU zeroes, except for 1 at the {k,£)th entry. 
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[i{eke + eek),z] = {zi - Zk)i{eki - e^) e P = ^n- (1.1.3) 

Incidentally, this implies that Cki + eik and i{eki + e^k) are two-dimensional eigenspaces^ 
of [adzY with eigenvalue {zi — ZkY- From (1.1.2) and (1.1.3) it follows that 

[A, z] = {zi - Zk) ^ {aki{eki + eik) + ihki{eki - eik)) ^ P = T-in (1-1-4) 

l<k<l<n 

and thus, for small A, we haveQ 

dM = d{e^ze-^) 

= d{z+[A,z] + . . .) 

n 

= Yi^^i n (^((^^ - ^fc)«fc^)'^((^^ - ^/c)&m), using (1.1.4) and (1.1.1) 



1 l<k<e<n 
n 



Y\_dziAl{z) Y\_ dakdbke- (1.1.5) 



l<k<l<n 



Therefore A^(z) is also the Jacobian determinant of the map M {z,U), such that 
M = UzU~^ G Tin, and thus dM admits the decomposition in polar coordinates: 

dM = Al{z)dzi . . . dZndU, U G SU{n). (1.1.6) 

In random matrix theory. Tin is endowed with the following probability, 

P(M G dM) = c„e-*"^(*^)dM, p{dz) = e~^^'^dz, (1.1.7) 

where dM is Haar measure (1.1.6) on 7i„ and c„ is the normalizing factor. Since dM 
as in (1.1.6) contains dU and since the probability measure (1.1.7) only depends on the 
trace of V{M), dU completely integrates out. Given E cM., define 

HniE) := {M G 7^„with all spectral points E E cR} cHn- (1.1.8) 

Then 

P(M G = [ c„e-^^WrfM = {^"t'ilSlfi 

Jhae) Jr" ^ (^) 111 P{dzk) 

As the reader can find out from the excellent book by Mehta it is well known 
that, if the probability P{M G dM) satisfies the following two requirements: (i) in- 
variance under conjugation by unitary transformations M i— UMU~^, (ii) the random 
variables Ma, 3?Mjj, SMj^, 1 < i < j < n are independent, then V{z) is quadratic 
(Gaussian ensemble). 



^adx{y) := [x,y]. 

'*A„(z) = Y\ i^i ~ ^j) is the Vandermonde determinant. 

l<z<j<n 



(ii) Symmetric ensemble 

Here we consider tlie non-compact symmetric spac^ SL{n,M.) / SO{n) witli a{g) = g^^^ 
. Then 

SL{n,R)/SO{n) = {gg'^ \ g G SL{n,R)} 

= {positive definite matrices with det = 1} 

with 

K = {ge SL{n, M) | aig) = g} = {g e SL{n, M) = g''} = SO{n). 

Then cr*(a) = — and the tangent space to G/K is then given by the space p = iS„ of 
symmetric matrices, appearing in the decomposition of s/(n, M), 

sl{n, R) = tQ)p = so{n) © iS„, i.e., a = ai + a2, ai G so{n), 02 G iS„ 
with Haar measure dM = Y\ dMij on iS„. 

l<i<j<ri 

A maximal abelian subalgebra a C p = iS„ is given by real traceless diagonal matrices 
z = diag(zi, ... , Zn)- Each M G p = 5„ conjugates to a diagonal matrix z 

M = e^z e"^, e K = SO{n), Ae so{n). 

A calculation, analogous to example (i) (1.1.5) leads to 

dM = I A„ (2) I rf^i . . . dzndU, U e SO (n) . 

Random matrix theory deals with the following probability on iS„: 

P{M G dM) = c„e-*"^(*^)c/M, p{dz) = e'^^'Uz, (1.1.10) 

with normalizing factor c„. Setting as in (1.1.8): Sn{E) C Sn is the subset of matrices 
with spectrum G E. Then 

P{M G 5„(E)) = / c„e-^^^(^^)rfM = {^"'t^i'S^fft (1.1.11) 

As in the Hermitian case, P{M G dM) is Gaussian, if P{M G dM) satisfies 

(i) invariance under conjugation by orthogonal conjugation M — > 0M0~^, (ii) Ma, Mij 

{i < j) are independent random variables. 

^Tlie compact version is given by SU{n)/ SO{n). 
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(iii) Symplectic ensemble 

Consider the non-compact symmetric spac^ SU*{2n)/USp{n) with a{g) 
where J is the 2n x 2n matrix: 



1 7-1 



and 



J :-- 



1 
-1 








1 
-1 








1 
-1 



V 



with J2 = -/, 



:i.i.i2) 



G = SU*{2n) = {ge SL{2n, C) \ g = JgJ'^}, 

K = {ge SU*{2n) \ a{g) = g} := Sp{n, C) n U{2n) 

= {ge SL{2n, C) | g^Jg = J} n {g e SL{2n, C) | g-' = f} 

= {geSL{2nX)\g-'=f^ndg = JgJ-'} 

= : USpin). 



Then, a* (a) = —JaJ J ^ and 



6 = {a G su*{2n) | a^,{a) = a] = sp{n, C) fl u{2n) 
= {a e C'"^'" I = -a, a = JaJ''} 



p = {a e su*{2n) | cr^(a) = -a) = su*{2n) fl iu{2n) 
= {ae C'"^'" I = a, a = JaJ^^} 



M= (MH)l<M<n,M; 



M 



Mi? M« 



with Ma = Mj^eC2^2 



= {self-dual n x n Hermitean matrices, with quaternionic entries} 

= '■ T2n- 

The condition on the 2x2 matrices M^i implies that M^k = Mkl, with Mk G M and 
the 2x2 identity I. Notice U Spin) acts naturally by conjugation on the tangent space 

^The corresponding compact symmetric space is SU{2n)/ Sp{n). 
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p to G/K. Haar measure on is given by 

n 

dM = \[dM, n dMifdMifdMji^dMi]\ (1.1.13) 

1 l<k<e<n 

since these Mij are the only free variables in the matrix M G A maximal abelian 
subalgebra in p is given by real diagonal matrices of the form z = diag(2;i, zi, Z2, Z2, Zn, Zn). 
Each M e p = can be written as 

M = e^ze"^, e K = USp{n), (1.1.14) 

with [ttke, bke, Cke, dke G M) 

^ = E ^'^^(4? - 4?) + &M(ea) + e«) + c,,(ei? - ej?) + 4.(el? + eg^) G ! 

l<k<t<n (1.1.15) 

in terms of the four 2x2 matrices [] 

,(0) _ ^ 1 M .(1) ^ \ .(2) - ^ M .(3) ^ ^ 

Since 



" = ' 1 ' " = ' " = -1 ' " V . 



■JO) 

^kl 






= (Zi - 


Zk){e^k} 




e P 




^ke 






= (Zi - 


Zk){e'2 




e P 




^ke 




A 


= {zi - 


Zk){,e^kl 


+ 4?) 


e p 






^ ^ik 5 


A 


= {Zi - 






e p, 


(1.1.16) 



[ 

[A, 2;] G p has the following form: it has 2x2 zero blocks along the diagonal and from 
(1.1.16) and (1.1.15), 

{{k, block in [A, z]) = {z, - Zk) ( ^ ^ ],ik<£). (1.1.17) 



-Cke + idke ake — ibkt 

Therefore, using (1.1.17), Haar measure dM on equals 

dM = d{e^ze-^) 

= d{I + A + ...)z{I~A + ...) 

= d{z+[A,z] + ...) 

= n '^^^ n d{{zi - Zk){aki + ihki))d{{zi - Zk){aki - ihki)) 

l<k<n l<k<e<n 

d{{zi - Zk){cki + idki))d{{zi - Zk){-Cke, + idkp)) 
= A'^{z)dzi- ■ ■ dzn Y\. ^dakdhtdckdd, 



IkE- 

l<k<e<n 



^e[,£ in (1.1.15) refers to putting the 2x2 matrix e^*^ at place (fc,£). 
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As before, define T2n{E) C as the subset of matrices with spectrum G E and 
define the probabihty: 

Remark: Notice is called the symplectic ensemble, although the matrices in p = 
are not at all symplectic; but rather the matrices in 6 are. 



1.2 Infinite Hermitian matrix ensembles 

Consider now he limit of the probability 

P{M G HniE)) = r A2f ^rr" m ^ ' ^^^"^ n / cg. (1.2.1) 

Dyson (see also Mehta [^) used the following trick, to circumvent the problem 
of dealing with oo-fold integrals. Using the orthogonality of the monic orthogonal 
polynomials pk = Pk{z) for the weight p{dz) on M, and the L^-norms hk = J^pl{z)p{dz) 
of the pkS, one finds, using (det A)^ = det{AA^), 

n 

A\z)l[p{dz,) 
1 

» n 

= / det{pi-i{zj))i<i,j<ndet{pk-i{ze))i<k,e<nY\p{dzk) 
"^K" k=i 

n „ 

= ^ {-'^T^'^'W I Pnik)-l{Zk)P7T'{k)-liZk)p{dZk) 
7r,7r'Go-„ k=l 



n— 1 „ n— 1 



in / p^(^)p(^^^)=n!n/^fc- (1-2.2) 



n 





For the integral over an arbitrary subset C M, one stops at the second equality, since 
the pnS are not necessarily orthogonal over E. This leads to the probability (1.2.1), 



P(M G HniE))) 



■ i-rn , — / det [ V pj_i{zk)pj-i{zg)] t\p{dzi) 

' l</c,^<n ^ 

1 /■ " 

— / dei{Kn{zk,zi))i<k/<n\\p{dzi), (1.2.3) 
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in terms of the kernel 

Pj-i{y) Pj-i{z 



K^{y,z):=J2'-^P^. (1.2.4) 



The orthonormahty relations of the Pk{y) / y/hk lead to the reproducing property for the 
kernel Kn{y,z): 

Kn{y,z)Kn{z,u)p{dz) = Kn{y,u), Kn{z, z)p{dz) = n. (1.2.5) 

Jr Jm. 

Upon replacing E"- by Yli dzi x in (1.2.3), upon integrating out all the remaining 
variables Zk+i, ■■■,Zn and using the reproducing property (1.2.5), one finds the n-point 
correlation function 

P(one eigenvalue in each [zi, Zi + dzi], i = 1, k) 

k 

l[p{dzi). (1.2.6) 

1 

Finally, by Poincare's formula for the probability P(ui?j), the probability that no 
spectral point of M belongs to E is given by a Fredholm determinant 

P{Menn{E')) = det(/-AfCf) 

oo „ k 

= l + ^(-A)M det(K^{z,,z,)) J[p{dz,), 

k=i J zi<...<zk ^ ' i<*j<«; ^ 

for the kernel K^{y, z) = Kn{y, z)Ie{z). 

• Wigner's semi-circle law: For this ensemble (defined by a large class of p's, in partic- 
ular for the Gaussian ensemble) and for very large n, the density of eigenvalues tends 
to Wigner's semi-circle distribution on the interval [—y/2n, \p2/n\: 



-V2n — z^dz, \z\ < v^2n 



density of eigenvalues 



vr 



0, \z\ > V2^. 



• Bulk scaling limit From the formula above, it follows that the average number of 
eigenvalues per unit length near z = ("the bulk") is given by \/2n/Tr and thus the 
average distance between two consecutive eigenvalues is given by Tr/\/2n. Upon using 



this rescaling, one shows (|^, ^ ^ 



lim ''!— Kn I —7==, I = ^"^^^^^ — ^ (Sine kernel) 
n/^y^ \^/2n V2nJ Tr{x - y) 
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and 

P(exactly k eigenvalues G [0, a]) = — — — ( — 1 det(J — \KI[o^a\ 



with 

det(/ - XKIio^a]) = exp / ' ' dx, (1.2.7) 

Jo ^ 

where /(x, A) is a solution to the following differential equation, due to the pioneering 
work of Jimbo, Miwa, Mori, Sato (' = d/dx) 



{xff = A{xf - f)i-f" - xf + f) , with fix; A) - for x ^ 0. 

TT 

(Painleve V) (1.2.8) 

• Edge scaling limit: Near the edge \^2n of the Wigner semi-circle, the scaling is \/2n^/^ 
and thus the scaling is more subtle: (see ^ 0, |4^, [6^ ) 



y = V2^ + ^ , , (1.2.9) 
and so for the kernel Kn as in (1.2.4), with the p„'s being Hermite polynomials, 

\ I . n . 1) \ 

lim — Kn V2n + — , V2n + — = K(u, v), 

where 

/»oo /»oo 

K{u,v)= / A{x + u)A{x + v)dx, A{u) = / e*'"'"'''/^rfx. 

Jo J-oo 

Relating y and u by (1.2.9), the statistics of the largest eigenvalue for very large n is 
governed by the function, 

P{X....<y) = P(^2rii(^-l) , forn/oo, 

= det(/ — -ft'/(_oo,«]) = exp ^— y {a — u)g'^{a)da 

with g{x) a solution of 

g" = xg + 2g'^ 



^ (Painleve II) (1.2.10) 



The latter is essentially the asymptotics of the Airy function. In section 5, I shall 
derive, via Virasoro constraints, not only this result, due to Tracy- Widom p4|, but 



14 



\1, p. 15 



also a PDE for the probability that the eigenvalues belong to several intervals, due to 
Adler-Shiota-van Moerbeke 0, |l2l. 

• Hard edge scaling limit Consider the ensemble of n x n random matrices for the La- 
guerre probability distribution, thus corresponding to (1.1.9) with p{dz) = z'^^'^e~^^'^dz. 
One shows the density of eigenvalues near z = is given by 4n for very large n. At this 



edge, one computes for the kernel (1.2.4) with Laguerre polynomials p„ |52, 30]: 



where K^^\u^v) is the Bessel kernel^ with Bessel functions J^,, 



xJy{xu)Jy{xv)dx 



2{u - v) 



Then 



P{ no eigenvalues G [0, x]) = exp J 



fin] 



du 



u 



with / satisfying 

{xf'f - A{xf - f)f + ((x - u^)f ~f)f = 0. (Painleve V) 



;i.2.11) 



;i.2.12) 



1.2.13) 



This result due to Tracy- Widom and a more general statement, due to ||TT] , will 
be shown using Virasoro constraints in section 5. 



1.3 Integrals over classical groups 

The integration on a compact semi-simple simply connected Lie group G is given by 
the formula (see Helgason 



f{M)dM 



G 



1 

\W\ .IT 



sm 



a{iH) 



dt I f{utu )du, t = e 
u 



H 



'1.3.11 



where A C G is a maximal subgroup, with q and a being the Lie algebras of G and A. 
Let du and dt be Haar measures on G, A respectively such that 



dt 



du = 1; 



A denotes the set of roots of q with respect to a; \W\ is the order of the Weyl group of 
G. 
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Integration formula (1.3.1) will he applied to integrals of f = e^i uttM^ ot;er the 
groups S0{2n), S0{2n+1) and Sp{n). Their Lie algebras (over C) are given respectively 
by 0„, b„, c„, with sets of roots: (e.g., see pO|) 



with 



A„, = {iecj, I <i <k, ±(ej + ej), ±(ej - e^), I <i < j <n}, 

£ = for c)„ = so{2n) 
£ = 1 for b„ = so{2n + 1) 
£ = 2 for c„ = sp{n). 

Setting H = i6, we have, in view of formula (1.3.1), 

a{iH) 



sm ■ 



dt 



Cn (ni<,<fc<„sin^sin^) UidOj for c)„ 

c.(ni,.<.,„sin^sin^-^^^'-— ^. 



2 , Hi sin ^dOj for b„, c„ 
dOj for dn 



c'^ n {cosOj - cos6'fc)^ < 

l<j<k'<n 



l<j<" 

n 

l<j<n 



1 — COS 9j 



dOj for bn 



JJ^ (1 — cos^ 6 j)d6j for 

, l<j<n 

for D„ 



l<i<n a/ 1 - 



4A2(z) n (1 



z, i — for 

^= for Cr, 



l-.l 



a = /? = -1/2 



for l)r 



c'^A2(^) n (1 - ^j)"(l + ^jO'^^'^j ^ a = 1/2, /3 = -1/2 for b. 



l<j<n 



a = /5 = 1/2 



for Cr 



For M G S'0(2n), Sp{n), the eigenvalues are given by e*^-' and e 1 < j < 
therefore, setting / = exp(^tfctrM'^) in formula (1.3.1), leads to 



(1.3.2) 
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where T„(z) are the Tchebychev polynomials, defined by Tn{cos6) := cosnO; in partic- 
ular Ti{z) = z. 

For M G S0{2n + 1), the eigenvalues are given by 1, e*^^ and e"*^^, 1 < j < n, 
which is responsible for the extra-exponential e^*' appearing in (1.3.2). 
Before listing various integrals, define the Jacobi weight 

Paf3{z)dz := (1 - z)"(l + zfdz, (1.3.3) 

and the formal sum 



1 

The arguments above lead to the following integrals, originally due to H. Weyl |]75|, and 
in its present form, due to Johansson [Q; besides the integrals over SO{k) = 0+(/c), 



the integrals over 0-{k) and U{n) will also be of interest in the theory of random 
permutations: 

f e5:rM.M^^^ = f A^{zrf[e^^^^^p,,_u{z,)dz, 

Jo{2n)+ J [-1,1]^ fj[ ^' ' 

f e^TutrM^dM = e^^*' ! A^izfUe^^'^^ pa _u{zu)dzk 

Jo(2n+lW 111 ^' ^ 



f e^rtarM^dM = [ A^{zrfle^^^''^pn,Jz,)dz,. 

« « 71 — 1 

/ e^^'^'^'^'dM = e^^^*,, / A^^,{zfY\e^^^^^p,..){zk)dz, 

[ e'^TutrM^dM = e^^^-^y*' [ A^{zff\ e^^'^^ p.i i.{z,)dz. 



f ^ETHUAP-s.mdM = - I |A„(^)|2freS?°(*>^^-^'0^ 
Ju{n) n\ J^s^)n 2mzk 

(1.3.4) 

1.4 Permutations and integrals over groups 

Let Sn be the group of permutations vr and S'2„ the subset of fixed-point free involutions 
7r°(i.e., (vr°)^ = / and TT^{k) ^ k for 1 < k < 2n ). Put the uniform distribution on 5"^ 
and ; i-e., all permutations or involutions have equal probability: 

P(vr„) = 1/n! and P«) = (1.4.1) 
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7r„ refers to a permutation in S'„, and 7r2,„ to an involution in S'^^. 

An increasing subsequence of n E Sn or is a sequence 1 < ji < ... < jfc < n, such 
that vr(ji) < ... < T^{ik)- Define 

L{T^n) = length of the longest increasing subsequence of 7r„ . (1-4.2) 

Example: for vr = (3, 1,4, 2, 6, 7, 5), we have L^ttj) = 4. 

Around 1960 and based on Monte-Carlo methods, Ulam conjectured that 

lim ^^^^ = Q exists. 



'n 

An argument of Erdos & Szekeres dating back from 1935 showed that E{Ln) > 
^a/^ ~ 1; and thus c > 1/2. In '72, Hammersley [^] showed rigorously that the limit 



exists. Logan and Shepp showed the limit c > 2, and finally Vershik and Kerov 



74| that = 2. In 1990, I. Gessel [31 showed that the following generating function is 



the determinant of a Toeplitz matrix: 

"2^ 



^ j-n / p'lTT \ 

-P{Ln < £) = det / e^^t^ose^i{k-m)6^Q _ ^^ 4 3) 

n=0 ^' ^"^0 / 0<k,m<e-l 



The next major contribution was due to Johansson [^| and Baik-Deift- Johansson [T7 



who prove that for arbitrary x G M, we have a " law of large numbers" and a " central 
limit theorem" , where F{x) is the statistics (1.2.10), 



lim — ^ = 1, and P ( — — -j^ < x] — > Fix), for n — > oo. 

n^oo 2^/n \ n^/° / 

A next set of ideas is due to Diaconis & Shashahani [^, Rains p7| , p8[ , Baik & Rains 
T8| . For a nice state-of-the-art account, see Aldous & Diaconis [|1^. An illustration 



is contained in the following proposition; the first statement is essentially Gessel's and 



the next statement is due to |26, 58, 18 



Proposition 1.1 The following holds 

°° j-n r _ 

V ^P(L(7r„) < £) = / e^^(^+^)dM (1.4.4) 



{/(£) 



[0,27r]' 



^<j<k<t l<k<i 



J2 i^P(L(7r°J <i)= [ e*^^dM. (1.4.5) 



n 

n=0 
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The proof of this statement will be sketched later. The connection with integrable 
systems goes via the following chain of ideas: 

Combinatorics 
i 

Robinson-Schensted-Knuth correspondence 
i 

Theory of symmetric polynomials 
i 

Integrals over classical groups 



Integrable systems 

All the arrows, but the last one, will be explained in this section; the last arrow will 
be discussed in sections 7 and 8. We briefly sketch a few of the basic well known facts 
going into these arguments. They can be found in MacDonald ||4^, Knuth Aldous- 



Diaconis [13] . Useful facts on symmetric functions, applicable to integrable theory, can 



be found in the appendix to Let me mention a few of these facts: 

• A Young diagram A is a finite sequence of non- increasing, non-negative integers 
Ai > A2 > ... > A^ > 0; also called a partition of n = |A| := Ai + ... + A^, with 
|A| being the weight. It can be represented by a diagram, having Ai boxes in the 
first row, A2 boxes in the second row, etc., all aligned to the left. A dual Young 
diagram A = (Ai > A2 > ...) is the diagram obtained by flipping the diagram A 
about its diagonal. 

• A Young tableau of shape A is an array of positive integers aij (at place {i,j) in 
the Young diagram) placed in the Young diagram A, which are non-decreasing 
from left to right and strictly increasing from top to bottom. 

• A standard Young tableau of shape A is an array of integers 1, ...,n placed in the 
Young diagram, which are strictly increasing from left to right and from top to 
bottom. The number of Young tableaux of a given shape A = (Ai > ... > Am) is 
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given by a number of formulae (for the Schur polynomial sx, see below)^ 

= ^{standard tableaux of shape A} 
= coefficient of X1X2 . . . x„ in sx{x) 

= TT , A = l^l'detf- -) (1.4.6) 

m ^ 

= l-^l' W i^i~^j)Wl~li with /ij := Ai — i + m, m := Ai. 

l<:i<j<m 1 

The Schur polynomial sx associated with a Young diagram A is a symmetric 
function in the variables xi,X2, ... (finite or infinite), defined by 

Sa(xi,X2,...) := Yl Yl^^^J■ (1-4-7) 

{(lij} tableaux of A ij 

The linear space A„ of symmetric polynomials in zi, x„ with rational coefficients 
comes equipped with the inner product 



If " 

-^^•^1^ l<k<l<n 1 



dzk 



2Tiizk 

f{M)g{M)dM. (1.4.8) 

U{n) 

An orthonormal basis of the space A„ is given by the Schur polynomials Sx{xi, a;„) 
in which the numbers aij are restricted to 1, n. Therefore, each symmetric func- 
tion admits a ^^Fourier series" 

f{xi,...,Xn) = ^ {f,Sx) Sx{Xi,...,Xn), with (sa, Sa') = ^aa' • (1.4.9) 

A with 

Ai<n 

In particular, one proves (see (1.4.6) for the definition of /^) 

(a;i + ... + x„)'== 5^/"sA, (1.4.10) 

jA| = fc 
\l<n 

If A = (Ai > ... > A£ > 0), withQ Xi = i > n, then obviously sa = 0. 



^h^j := Xi + Xj — i — j + 1 is the hook length of the i, jth box in the Young diagram; i.e., the length 
of the hook formed by drawing a horizontal line emanating from the center of the box to the right and 
a vertical line emanating from the center of the box to the bottom of the diagram. 

^Remember, from the definition of the dual Young diagram, that Ai = 
the length of the first column of A 
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Robinson-Schensted-Knuth correspondence: There is a 1-1 correspondence 



(P, Q), two standard Young 
tableaux from 1, . . . ,n, where 
P and Q have the same shape 



Given a permutation ii, ...,in, the correspondence constructs two standard Young 
tableaux P, Q having the same shape A. This construction is inductive. Namely, 
having obtained two equally shaped Young diagrams Pk,Qk from ii, ...,ik, with 
the numbers (ii, ifc) in the boxes of P^ and the numbers {1, ...,k) in the boxes 
of Qk, one creates a new diagram Qk+i, by putting the neoii number ik+i in the 
first row of P, according to the following rule: 

(i) if ijfc+i > all numbers appearing in the first row of P^, then one creates a new 

box with ife+i in that box to the right of the first column, 

(ii) if not, place ik+i in the box (of the first row) with the smallest number 
higher than ik+i- That number then gets pushed down to the second row of 
Pfc according to the rule (i) or (ii), as if the first row had been removed. 

The diagram Q is a bookkeeping device; namely, add a box (with the number 
+ 1 in it) to Qk exactly at the place, where the new box has been added to Pfe. 
This produces a new diagram Qk+i of same shape as Pk+i- 

The inverse of this map is constructed essentially by reversing the steps above. 



Example: n = (5, 1, 4, 3, 2) G 5*5, 



1 4 
5 



1 3 
2 



//I 2\ 

3 
4 

VVs J 

and so L^lir) = 2 = ^columns of P or Q. 



Hence tt — > (P(7r), g(7r)) 



1 3 
4 



1 3 
2 



/I 3\\ 

2 



4 



V5 J 



I 



1 2 
3 
4 
5 

1 3 
2 
4 
5 



The Robinson-Schensted-Knuth correspondence has the following properties 
• TT ^ (P, Q), then tt"^ ^ (Q, P) 
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• length (longest increasing subsequence of tt) = # (columns in P) 

• length (longest decreasing subsequence of vr) = 7^ (rows in P) 

• Tc'^ = I, then 71 ^ (P, P) 

• TT^ = /, with k fixed points, then P has exactly k columns of odd length. ^ ^^-^ 
From representation theory (see Weyl and especially Rains [0), one proves: 



Lemma 1.2 The following perpendicularity relations hold: 



(i) / sx{M)s^{M)dM ={sx,s^)= 6x^ 

Ju(n) 



(ii) / sx{M)dM = 1 for X = {Xi > . . . > Xk > 0), k < n, Xi even 

Jo{n) 

= otherwise. 



(Hi) / sx{M)dM = 1 for Xi even, Xi < 2n , 

J Spin) 

= otherwise. (1.4.12) 
Proof of Proposition 1.1 : On the one hand, 

((Xi + . . .+Xnf, (Xi + . . . + 

I A I — I /i| — fc 

A ]^ , /i ]^ < 71 

= E(/' 

|A|^fc 
A-j^ <n 



eX\2 



eX\2 



= #{(-P, Q), standard Young tableaux, each of arbitrary shape A 

with |A| = k, Ai < n} 
= ^{-Kf^ e Sk such that L^iTk) < n}. (1.4.13) 

On the other hand, notice that, upon setting 6j = 6j + 61 for 2 < j < n, the expres- 
sion Y\ ~ e*^*"!^ is independent of 9i. Then, setting Zk = e*^*, one computes: 

l<j<k<n 
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{(Xi + ... + Xn)'', (xi + ... + XnY) 

= 1 / e^^^Hl + 4 + ■ ■ ■ + z:fe-^''\l + 4 + ■ ■ ■ + 

l<j<fe<n 

upon setting 9j — 9j + 6i , for j > 2 and z'f. — e*^fc , 

1 /"^'^ 

= — / e^M^M^i X (an n - 1-fold integral) 
^! Jo 

= 6ke{{xi + . . . + Xn)Mxi + . . . + Xnf) = 6ke [ {TrMf'^dM. (1.4.14) 
It follows that 

{tr{M + M)fdM 
E f ■ ) / {{trMy{t^r^)dM 



0<jr<A; 



0, if k is odd (because then j ^ k — j for all < j < k) 

) ItrMl^dM, if A; is even. ^^"^-^^^ 

Combining the three identities (1.4.13), (1.4.14) and (1.4.15) leads to 

#{7rfc e Sk such that L(7rfe) < n} = ^ ^ (Tr(M + M)f^dM. (1.4.16) 

Finally 
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ni 

n=0 



n=0 



E^/ (tr(M + M)fW 



1 f pVt{zi+z^^+...+ze+z-^) TT _ iefc|2 TT ^ 



where = 6*^*= , 



1^ ,jefc|2 TTg2ytcos6li,^ 



showing (1.4.4) of Proposition 1.1. The latter also equals: 

dzk 



= 1 / A,(^)A,(^-)nfe^(^'=+^-'=) 



2m Zk 



dzk 

zz 'z-C -e-'-'-^'—r- 
" " 2mzi 



k-l-m~l Vi{zk+Zk) 



detl / (,'^Vico.e^Kk-m)e^Q 



l<k,m<e 



confirming Gessel's result (1.4.3). 

The proof of the second relation (1.4.5) of Proposition 1.1 is based on the following 
computation: 
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[ {TiMfdM = y^f^f sx{M)dM, using (1.4.10), 

-^Oin) i;,,^, Jo{n) 

<n 

, using Lemma 1.2, 



Xi<n 

Xi even 

— y2 J using duality, 

Ai<n 

Ai even 

_^ f {P,P),P standard Young tableau of shape A 
\ with |A| = fc, Ai < n, Aj even 

— # {tt^ e 5"°, no fixed points and -L(7r°) < n} . 

(1.4.17) 

In the last equality, we have used property (1.4.11): an involution has no fixed points 
iff all columns of P have even length. Since all columns Aj have even length, it follows 
that |A| = k is even and then only is /^^^^(Tr M)*^(iM > 0; otherwise — . Finally, one 
computes, 

fe=0 

°o f2k 2^h\ 

= E^(2^#{^2fee'52V ^«)<^}, using (1.4.1), 



fc=0 



oo 
fc=0 

oo ,^ ~ 

V— / (Tr M)^dM, using (1.4.17), 
^' Join) 



e''''-^dM, 

0(n) 



ending the proof of Proposition 1.1. 
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2 Integrals, vertex operators and Virasoro relations 

In section 1 , we discussed random matrix problems over different finite and infinite ma- 
trix ensembles, generating functions for the statistics of the length of longest increasing 
sequences in random permutations and involutions. One can also consider two Hermi- 
tian random matrix ensembles, coupled together. All those problems lead to matrix 
integrals or Fredholm determinants, which we list here: (/9 = 2, 1, 4) 

• / e-^^W^M = c„/ \Ar,{z)ff\pizk)dzk 

JHn{E), Sn{E) or Tu{E) J E^^ 1 

HUE1XE2) 

e^'^^dM 




L 



0{n) 
JUin) 

det(/ - XK{y, z)Ie{z)) with K{y, z) as in (1.2.4). ^2 1) 



The point is that each of these quantities admit a natural deformation, by inserting time 
variables ii, ^2, ••• and possibly a second set Si, S2, in a seemingly ad hoc way. Each 
of these integrals or Fredholm determinant is then a fixed point for a natural vertex 
operator, which generates a Virasoro-like algebra. These new integrals in ti,t2.... are 
all annihilated by the precise subalgebra of the Virasoro generators, which annihilates 
Tq. This will be the topic of this section. 



2.1 Beta-integrals 

2.1.1 Virasoro constraints for Beta-integrals 

Consider weights of the form p{z)dz := e^^'^^^dz on an interval F = [A, B] C M, with 
rational logarithmic derivative and subjected to the following boundary conditions: 

-^^V'^^ = hm / = for all k > 0, (2.1.1) 

and a disjoint union of intervals, 

2r 

E^[j [c2i-i,C2i]c F CR. (2.1.2) 
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These data define an algebra of differential operators 

O 



(2.1.3) 



Take the first type of integrals in the list (2.0.1) for general /3 > 0, thus general- 
izing the integrals, appearing in the probabilities (1.1. 9), (1.1. 11) and (1.1.18). Con- 
sider t-deformations of such integrals, for general (fixed) /? > 0: (t := (ti,t2, •••); 
c = (ci, C2, C2r) and Z = {Zi, Zn)) 



J„(t,c;/5) : = 



« n 

/ \A^{z)f 11 fe^^'^'l p{zk)dzk) for n>0. (2.1.4) 
■^E" k=i ^ ^ 

The main statement of this section is Theorem 2.1, whose proof will be outlined in the 
next subsection. The central charge (2.1.9) has already appeared in the work of Awata 



et al. ig. 

Theorem 2.1 (Adler-van Moerbeke 0, P) The multiple integrals 



/„(t,c;/3) : = 



„ n 

/ \/\^{z)\l'\\(e^^'^^^p{zk)dzu) , for n>0 
k=i ^ ^ 



(2.1.5) 



and 



4(t,c;-):= / \A^iz)\'/^l[(e^^'^^'^p{zk)dzX for n > 0, (2.1.6) 
with Iq = 1, satisfy respectively the following Virasoro constraint^ for all k > —1: 
-B, + (a, '^jfi.nli, n) - h ^tfi+i,n(i, n)) ) /„(t, c; /3) = 0, 



2 ' /3 ^ 2 2 ' /? 

(2.1.7) 

in terms of the coefficients ai, hi of the rational function (— log p)' and the end points Ci 
of the subset E, as in (2.1.1) to (2.1.2). For aline Z, the ^ji^i(t,n) and ^4^i(t,n) 
form a Virasoro and a Heisenberg algebra respectively, interacting as follows 

'k^ - k' 



/3¥(2) /3^{2) 
'^fc,n' '"£,11 



1) PjW 

k,n^ '"£,12 



(k-i) 



j(2) 



+ C 



12 



^k-e 



k 



P 



Sk,-e, 



(2.1.8) 



10 



When E equals the whole range F, then the Bk's are absent in the formulae (2.1.7). 
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with central charge 
Remark 1: The are defined as follows: 

"jS = f E ^ "Jl'i '^>l ^ + (i - f ) ((* + 1) - "S) ■ (2.1.10) 

i+j=k ^ ^ 

Componentwise, we have 

'^^nit. n) = ^4^) + njf and ^4°) = njf^ = n<^o. 

and 
where 

i+j=k ■' —i+j=k ■' —i—j=k 

We put n explicitly in ^lf'n(^,n) to indicate the nth component contains n explicitly, 
besides t. Note that for /3 = 2, (2.1.10) becomes particularly simple: 

/3|(2) _ . 2^(1) 2^(1) . 

Remark 2: The Heisenberg and Virasoro generators satisfy the following duality prop- 
erties: 



'jg(M.) = -f *j£(-f.-|).«>0. (2.1.12) 

In (2.1.7), '^J^^^ (— /3t/2, — 2n//3) means that the variable n, which appears in the nth 
component, gets replaced by 2n/f3 and t by —f3t/2. 

Remark 3: Theorem 2.1 states that the integrals (2.1.5) and (2.1.6) satisfy two sets of 
differential equations (2.1.7) respectively. Of course, the second integral also satisfies 
the first set of equations, with P replaced by 4//3. 
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2.1.2 Proof: /3-integrals as fixed points of vertex operators 

Theorem 2.1 can be established by using the invariance of the integral under the trans- 
formation Zi I— > Zi + ef{zi)z^^^ of the integration variables. However, the most transpar- 
ent way to prove Theorem 2.1 is via vector vertex operators, for which the /5-integrals 
are fixed points. This is a technique which has been used by us, already in 0. In- 
deed, define the (vector) vertex operator, for t = (ti,t2, •••) ^ u E C, and setting 
Xiz) := {l,z,z^,...). 



which acts on vectors f(t) = {fo(t), fi{t), ...) of functions, as followsQ 



(2.1.13) 



(X^(t,n)/(t)), 



(|nn"-Vn-i(t-/3[n-^]). 



For the sake of these arguments, it is convenient to introduce the following vector 
Virasoro generators: (t) := ( „(i, ^))nez- 

Proposition 2.2 The multiplication operator z^ and the differential operators ■§^z^^^ 
with z G C*, acting on the vertex operator Xjsit, z) , have realizations as commutators, 
in terms of the Heisenberg and Virasoro generators ^J^^''(t, ra) and ^I^^\t,n) : 



^J«(t),X,(t,z) 
^lf\t),X,{t,z) 



Corollary 2.3 Given a weight p{z)dz on R satisfying (2.1.1), we have 
^^z'^-''f{z)Xp{fz)p{z) = 

Li>0 



E - h 'tU,{t)),x,{fz)p{ 



(2.1.14) 



(2.1.15) 



Proof: Using (2.1.13) in the last line, compute 
d 



z^+^f{z)Xp{t,z)p{z) 



dz 



^/(^)) z^^^Xp{t,z)p{z)+p{z)^^ {z'^'fiz)X,it,z)) 



OO 



- I Y.kz'^^^'X,{t,z)j p(z) +p(z)^ I 

OO 

J2k ^3'il^„Xp{t,z)p{z 





aiZ' 



fc+i+i 



Xg{t, z) 



+ 



X^a, Plfl,Xp{fz)p{2 



(2.1.16) 



^^For a S C, define [a] := (a, ^,^,...) S C°°. The operator A is the shift matrix, with zeroes 
everywhere, except for I's just above the diagonal, i.e., {Kv)n = Vn+i- 
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establishing (2.1.14). ■ 

Giving the weight pE{u)du — p{u)lE{u)du, with p and E as before, define the 
integrated vector vertex operator 



^(3{t,PE) ■= / dup{u)Xp{t,u). 



(2.1.17) 



E 



and the vector operator 



2r 



1 * j>0 

consisting of a c-dependent boundary part Bk and a (t, n)-dependent Virasoro part Vk- 
Proposition 2.4 The following commutation relation holds: 



[Dk,Yp{t,pE)]^0 
Proof: Integrating both sides of (2.1.14) over E, one computes: 

ljzj-{z''^'f{z)Xf,{t,z)p{z)) = ^(-1)X'+V(q)XM^,cOp(q) 



(2.1.19) 



2r 



d_ 

^ dCiJE 

[Bk,Yp{t,pE)], 



while on the other hand 



dz 



E 



,i>0 



"tfi - b^ ^tfi+i) , / dzpE{z)Xp{t,z) 

= [Vfe,Y^(i,PE)]. 
Subtracting both expressions (2.1.19) and (2.1.20) yields 

= [Sfc - Vk, Yp{t, pe)] = [Vk, Yp{t, pe)] , 
concluding the proof of proposition 2.4. 



(2.1.20) 



(2.1.21) 
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Proposition 2.5 The column vector, 

V-^^" fe=l / n>0 

is a fixed point for the vertex operator Yp{t,pE): (see definition (2.1.17)) 

{Y^{t,pE)I)^^In, n>l. (2.1.22) 

Proof: We have 

» n 

= / \A^{z)fT\(e^r'^^llE{z,)p{^,)dz,) 

fc=i ik=i 
dupE{u)e^^'^^'\uf^''-'^ 

- „ n— 1 

g-/3Erv4 / \K-i{z)\^\\{e^^'^''hE{z,)dz,) 
/ dupE{u)\u\'^^''-^^e^^'''''e~^^^^^^Tn-l{t) 
[Yp{t,pE)I{t)) . (2.1.23) 



Proof of Theorem 2.1 : Prom proposition 2.4 it follows that for n > 1, 

= [Vk,(Yp{t,pE)r]i 

= VkY^it, peTI - Y^it, PeYVJ. (2.1.24) 

Taking the nth component for n > 1 and A; > — 1, and setting X^(t, u) — **"'e '^^~rmi 
and using (2.1.21) 

= {VkI-Yp{t,pErVkI)^ 

= {Vkl)n- j dupE(u)Xp(t;u)(\ufr-\..j dupE(u)X^(t;u)(VkI)o 
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Indeed {VkI)o = for A; > —1, since Tq — 1 and involves ^j'f'-, ^J^'' ^-^d jf' for 
k > -1: 

^jj^^ is pure differentiation for k > —1; 
^jj.^hs pure differentiation, except for k = —1; but 

appears with coefficient n(3, which vanishes for n = 0; 
jj!^^ appears with coefficient n((n — 1)| + 1), vanishing for n = 0. 



2.1.3 Examples 

Example 1 : Gaussian /3-integrals 

The weight and the Oj and 6j, as in (2.1.1), are given by 

Qq = l,bo = 0, bi = 2, and all other Oj, hi — 0. 
From (2.1.5), the integrals 

/„= / \A^{z)ff[e-^'^+^r=^'^^idz, (2.1.25) 
fe=i 

satisfy the Virasoro constraints, for A; > —1, 

-B.In = - E = (- 'jSl.n - « + ^tfi.n) 4- (2-1.26) 

Introducing the following notation 

<Ji = {n-'-^)l3 + i + l-bo={n-'-^)p + i + l, (2.1.27) 
the first three constraints have the following form, upon setting F„ = log/„. 

For later use, take the linear combinations 

Pi = -i^_i, V2 = -^Bo, 1)3 = (^1 + yS-i) , (2.1.29) 
such that each Dj contains the pure term dF/dti, i.e., "DjF = + ... 
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Example 2 (Laguerre /3-integrals) 

Here, the weight and the a, and 6j, as in (2.1.1), arc given by 

f z 

oq = 0, 0,1 = 1, 6o = —O'^ bi = 1, and all other a^, bi = 0. 
Thus from (2.1.4), the integrals 



» n 

4=/ \A^{z)fl[zte-^^+^^=Mdz, (2.1.30) 
fc=i 



satisfy the Virasoro constraints, for k > —1, 



2r r. 

-BJ^ = -J2 c^^/n = (- ^Jfi, - a + ^J«,„) 4. (2.1.31) 

Introducing the following notation 

Gi = =(n-^^)/3 + i + l-6o = (ri-^^)/3 + i + l + a. 



the first three have the form, upon setting F — — log/„. 

Again, replace the operators Bi by linear combinations "Dj, such that PjF = |^ + 

= ©2 = --Bo - (7iB_i, ©3 = --Bi - (72^0 - (7iC72-B_i. (2.1.32) 

Example 3 (Jacobi /5-integral) 

This case is particularly important, because it covers not only the first integral, but 
also the second integral in the list (2.0.1), used in the problem of random permutations. 
The weight and the and bi, as in (2.1.1), are given by 

p(z) := e-" = (1 - zm + zf, V" = f = "-"+'7"'' 
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ao — l,ai — 0, a2 — —1, bg — a — b,bi — a + b, and all other a^, bi — 0. 
The integrals 

p n 

/ \/^n{z)fT\{l- ZkYil^ zu)'e^r=.tAdzj, (2.1.33) 
satisfy the Virasoro constraints {k > —1): 

1 * 
= ( - + ^0 ^tfi, + 6i 4. (2.1.34) 



Introducing cTj = (n — ^)/3 + ^ + 1 + 6i, the first four have the form (k — —1, 0, 1, 2) 



9*3 dt2 dti ^ 9ii+3 dti+i dtidt2 

dF OF 

' oti dta dt2 ^ dti+4 dti+2 



/3 9^ 



(3 f OF ^ _ dF ^ dFdF\ 



(2.1.35) 



2.2 Double matrix integrals 

Consider now weights of the form 

p{x,y) = e^'.^>^'^^^'y'p{x)p{y), (2.2.1) 
defined on a product of intervals Fi x F2 C R^, with rational logarithmic derivative 



i>0 '^iV 
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satisfying 

lim f{x)p{x)x'' = lim f{y)p{y)y^ = for all k>0. (2.2.2) 

x->dF y~>dF 

Consider subsets of the form 

E = EixE2:= U[=i[c2.-i, C2^] X UU[c2i-i, C2i] C Fi X F2 C (2.2.3) 

A natural deformation of the second integral in the list (2.0.1) is given by the following 
integrals: 



hit, s,r;E) = [ [ A„(x) A„(t/) TT (^Zrit^-l-^^vl 
J k=i 



'^^''^p{xk,yk)dxkdyk (2.2.4) 



In the theorem below, J^?^^ and are vectors of operators, whose components are 
given by the operators (2.1.10) for /5 = 1; i.e.. 



/3=1, t^-s 



thus, from (2.1.10) and (2.1.11), one finds: 

JS(^) = I {4"\t) + {2n + k + 1) J«(t) + n{n + 1) ) , (2.2.5) 

satisfying the Heisenberg and Virasoro relations (2.1.8), with central charge c = —2 and 
c' = 1/2. 

The tti, di, bi, bi, Ci, Ci, Vij given by (2.2.1), (2.2.2) and (2.2.3) define differential oper- 
ators: 



1 * i>0 m,e>l m+k+i/ 



i>0 m,e>l --mZ+k+i 



(2.2.6) 

Theorem 2.6 (Adler-van Moerbeke 0, §]) Given p{x,y) as in (2.2.1), the integrals 

In{t, s, r; E) := [[ A„(x) A„(t/) TT e>^^^^'^^l-'^y^^ pixk, yk)dxudyk (2.2.7) 

k=i 

satisfy two families of Virasoro equations for k > —1: 

Vk,nIn{t,s,r;E) = and AMt, s,r; E) = 0. (2.2.8) 
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The proof of this statement is very similar to the one for /5-integrals. Namely, define 
the vector vertex operator, 



_d vZlJL 



X{uv) 



(2.2.9) 



which, as a consequence of Proposition 2.2, interacts with the operators J^^t) 



as follows: 



'u'^Xi2(t, s; M, v) 

d 



4')(t),Xi2(t,s;M,t;) 
jf (t),Xi2(t,s;M,t;) 



(2.2.10) 



A similar statement can be made, upon replacing the operators and ^^i*^^^ by 

and -^v^^^^ and upon using the J^*^(s)'s. 

Finally, one checks that the integral vertex operator 

Y{t,s;pE) ■■= J j dxdyp{x,y)Xi2{t,s;x,y) (2.2.11) 

commutes with the two vectors of differential operators = {T>k,n)nez, as in (2.2.6): 



Vk,Y{t,s-pE) = Vk,Y{t,s;pE) 



0, 



and that the vector / = (Jq = 1, h, ...) of integrals (2.2.7) is a fixed point for Y(t, s; pe), 

Y{t,s;pE)Iit,s,r;E) = I{t,s,r;E). 

Then, as before, the proof of Theorem 2.6 hinges ultimately on the fact that T>kfi 
annihilates /q = 1. 

2.3 Integrals over the unit circle 

We now deal with the fourth type of integral in the list (2.0.1), which we deform, this 
time, by inserting two sequences of times ti, ^2, ••• and Si, S2, .... The following theorem 
holds: 

Theorem 2.7 (Adler-van Moerbeke [0) The multiple integrals over the unit circle , 

dzk 



(51)" 



\A^{z)\'lle^^^'^^l-^^^^'^^, n>0, 



k=l 



27rizk ' 



(2.3.1) 
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with Iq — 1, satisfy an SL(2,'L)-algehra of Virasoro constraints: 

9 = 

"D^ ^J„(t, s) = 0, for ^ k = 0, 9 arbitrary ^ only, 

k = l, 9 = 1 



(2.3.2) 



where the operators Vf, ,^ := Vl^^{t, s,n), A; e Z, n > are given by 

Kn ^klit, n) - J2,„(-s, n) - k n) + (1 - 9)/X^{-s, n)) , (2.3.3) 

withStlit,n) := ^S^lit,n)\^^^, as m (2.1.11). 
The explicit expressions are 

V-^In = K](i + l)ti+i^-^(i-l)si_i^+n(ti + ^ ) I /„ = 



, i>l 



i>2 



dsi 



vjr. = (|-» + i)^m|: + E(i-i)i.-i|: + n(s, + A)j 



(2.3.4) 



/„ = 0. 



Here the key vertex operator is a reduction of X^it, s;u,v), defined in the pre- 
vious section (formula (2.2.9)). For all k E Z,, the vector of operators Vf.{t,s) 
(T>1 j^{t, s,n)^ form a realization of the first order differential operators 
ing the vertex operator Xi2(t, s; u, u~^), namely 



du 



U ' % US- 



du u 



Xi2{t,s;u,u ) 



u 



(2.3.5) 



Indeed, 

u—u^'Ki2{t, s; u, u~^) 
du 

= ( ^u^+^ - ^v^-'' - k9u'' - k(l - 9)v-''] Xi2(t, s; u, v) 
\du dv J 

= (t) - jS(-s) - k (9I^^\t) + (1 - 9)l'i\-s)) ,X,2(t,s;u, -u)] 

= [Vl{t,s),X,2{t,s;u,u-']. 

The Vf, := V^^t, s) satisfy Virasoro relations with central charge = 0, 

[VIV',] = {k - i)vU„ (2.3.6) 
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and, from (2.3.5) the following commutation relation holds: 

[Vl{t,s),Y{t,s)]=0, with Y(t,s) := / -^X,2{t, s;u,u-'). (2.3.7) 

Jgi 2mu 

The point is that the column vector I{t, s) = (Jq, h, ...) of integrals (2.3.1), is a fixed 
point for Y{t, s): 

(Y(t,s)J)„ = /„, n> 1, (2.3.8) 
which is shown in a way, similar to Proposition 2.5. 

Proof of Theorem 2.7: Here again the proof is similar to the one of Theorem 2.1. 
Taking the nth component and the nth power of Y(t, s), with n > 1, and noticing from 
the exphcit formulae (2.3.4) that (Vl{t, s)!)^ = 0, we have, by means of a calculation 
similar to the proof of Theorem 2.1, 

= ([I?^,Y(t, .)"]/)„ 

= {viY{t,sri-Y{t,srvti)^^ 

= {Vll-Y{t,srvll)^={Vll)^. ■ 



3 Integrable systems and associated matrix inte- 
grals 

3.1 Toda lattice and Hermitian matrix integrals 

3.1.1 Toda lattice, factorization of symmetric matrices and orthogonal poly- 
nomials 

Given a weight p{z) = e"^*-^-* defined as in (2.1.1), the inner-product over C M, 

{f,9)t= [ f{z)g{z)pt{z)dz, withpi(z):=eS*'^V(^), (3.1.1) 



leads to a moment matrix 

mn{t) = {fiij(t))o<ij<n = {{z\z^)t)o<i,j<n, (3.1.2) 

which is a Hdnkel matriJQ, thus symmetric. This is tantamount to Anioo = niooA^, 
where A denotes the shift matrix; see footnote 11 . As easily seen, the semi- infinite 
moment matrix m^o evolves in t according to the equations 

dfiij , ^_ dm^ ^ / commuting \ /o i on 

- pi+kj, and thus — — = A m^. ^ , (3.1.3) 



^ „.oo. , ^g^^^^ ggj^g 



Hankel means /lij depends on i + j only 
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Another important ingredient is the factorization of moo into a lower- times an upper- 
triangular matrix^ 

fnoo{t) = S{t)~^S{t)~^~^, S{t) = lower triangular with 

non-zero diagonal elements. 

The main ideas of the following theorem can be found in 0, Remember c = 
(ci, C2r) denotes the boundary points of the set E. dM refers to properly normalized 
Haar measure on Tin- 

Theorem 3.1 The determinants of the moment matrices 

r„(t, c) := detm„(t,c) = — Al{z)Y\_Pt{zk)dzk 



k=l 



^tri-ViM)+ETUM^)dM, 



Hn{E) 



satisfy 

(i) Virasoro constraints (2.1. 7) for (3 = 2, 



1 * i>0 



ai J 



(2) 



(3.1.4) 



)jr„(t,c) = 0. (3.1.5) 



(ii) The KP-hierarch-^ (k = 0, 1, 2, ... j 

, d I d I d . 



1 



of which the first equation reads: 



d_ 



+ 3 



_d_ 



-4- 



dtidts 



logr„ + 6 



2 dtidtk+3 



rnOTn = 0, 



0. 



(3.1.6) 



(iii) The standard Toda lattice , i.e., the symmetric tridiagonal matrix 

/ « . _ /__ \V2 



L{t) := S{t)AS{t) 



-1 



TQT2 



1/2 







i-^^^f. (^) 

log 



1/2 



T1T"3 



1/2 



dtl 



\ 



(3.1.7) 



^•^This factorization is possible for those t's for which r„(t) := detm„(t) ^ for all n > 0. 
^^Given a polynomial p(ti, t2, ...), define the customary Hirota symbol p{dt)f o g 
p{-^, ...)f(t + y)g(t — y) . The s^'s are the elementary Schur polynomials e^i°*'^' 

E^>oS*(0^* and Si{d) S£( J^, i J^, . . . ). 
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satisfies the commuting equation^ 

dL_ _ 
dtk 



(3.1.8) 



(iv) Eigenvectors of L: The tridiagonal matrix L admits two independent eigenvectors: 

• p{t;z) = {pn{t] z))n>o satisfying {L{t)p{t; z))n = zpn{t;z), n > 0, where 
Pn{t', z) are nth degree polynomials in z, depending on t E C°°, orthonormal 
with respect to the t-dependent inner produc^\ (3.1.1) 

{pk{t; z),pi{t; z))t = Skf, 

they are eigenvectors of L, i.e., L{t)p{t; z) = zp{t; z), and enjoy the following 
representations: = iX^ z"^ ^ ...Y) 



Pnit;z) 



iSit)x{z))n 
1 

\/rn{t)Tn+i{t) 



( 



det 



m„(t) 



z'^K 



\ fJ'nfl ■ ■ ■ f^n,n-l 



1 \ 

Z 



-l/2 'rn{t-[z ^]) ^ 



(3.1.9) 



q{t,z) = {qn{t;z))n>o, with quit] z) := z ^^^^pt{u)du, satisfying {L{t)q{t; z 
zqn{t; z), n > 1; qn{t; z) enjoys the following representations: 



qn{t;z) := z 



n Z — U 



pt{u)du = {S^-\t)x{z-'))^ 

= {S{t)mUt)xiz-'))^ 

Tnit) ■ 

(3.1.10) 

In the case (3 = 2, the Virasoro generators (2.1.11) take on a particularly elegant 
form, namely for n > 0, 



= E^JSWJSW^ = Ji'\t) + 2nji'\t)+n'6^ 

i+j=k 

4!iW = Jl'\t)+n5. 



Ok 



JOk, 



means: take the skew-symmetric part of (); for more details, see subsection 3.1.2. 
^^The explicit dependence on the boundary points c will be omitted in this point (iv). 
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withQ 

i-\-j=k ■' —i+j=k •' —i—j=k 

Statement (i) is already contained in section 2, whereas statement (ii) will be established 
in subsection 3.1.2, using elementary methods. 



3.1.2 Sketch of Proof 

Orthogonal polynomials and r-function representation: The representation 
(3.1.4) of the determinants of moment matrices as integrals follows immediately from 
the fact that the square of a Vandermonde determinant can be represented as a sum of 
determinants 

A2K,...,«„)= 5^det 



l<k,l<n 



Indeed, 



n!r„(t) = n!detm„(t) 

^ det ( / z^t^-^ pt{z^(^k))dz^{k) 



l<k,i<n 



J2 f det 

Je" \ ^ ' / l<k,e<n 

n 

^l{z)^Pt{zk)dzk, 



^" k=i 



whereas the representation (3.1.4) in terms of integrals over Hermitian matrices follows 
from section 1.1. 

The Borel factorization of moo is responsible for the orthonormality of the polyno- 
mials pn{t; z) = {S{t)x{z))n; indeed, 

{Pk{t; z),pe{t; 2;))o<fc/<oo = / Sx{z){Sx{z)y pt{z)dz = SrUooS'^ = I. 

J E 

Note that Sx{,z){Sx{,z)Y should be viewed as a semi-infinite matrix obtained by multi- 
plying a semi-infinite column and row. The determinantal representation (3.1.9) follows 
at once from noticing that {pn{t; z) , z^) = for < /c < n — 1, because taking that 
inner-product produces two identical columns in the matrix thus obtained. From the 
same representation (3.1.9), one haspn(t; z) = h^^z"' + ... , where hn := (Tn+i/T„(t))^''^. 

^'''The expression jj^^^ = for k = 0. 



41 



§3, p. 42 



The "Sato" representation (3.1.9) of z) in terms of the determinant Tn{t) of the 
moment matrix can be shown by first proving the Heine representation of the orthogonal 
polynomials, which goes as follows: 

hnPn{t;z) 



( 



1 , 

— det 



m„{t) 



1 \ 

z 







\ l^n,0 






^ / 














/ 


n 


1 


rt 1 

11 


1 


\ 


\ 










u\ 


ul 


"n 


2; 




n 


f 


det 














W^pt{ui)dui 














,,2n-2 
"n 


2" 


-1 


1 








V 






.,2n-l 
"n 






/ 










ul 


ul ■■■ 


"n 


1 


\ 














ul 




z 








f 


det 














ulul ■■■ul~^ 


T„ 


J En 








■■■ 






1 










I 




u^ ■■■ 


< 




/ 










/ 


<m 


^l{2) ■ • ■ 

K{2) ■ ■ • 




1 


\ 












Z 






f 


det 
















Tn 


Je"^ 








7/"^^ • • • 






1 










V 




K{2) ■ • ■ 


K{n) 




/ 





7/0 7/^ 
^a(l)^o-(2) 



■ ■ M 



1 

'^n Je" 



det 



ul 
u\ 


ul ■ 

u\ 


■ < 

■ < 


z 


u'l-^ 


uT^ ■ 






ul 


ul ■ 


■ K 





1 

for any permutation a & Sn 
\ 



1 

If " 

- — / A^(u) \ \{z — Uk)pt{uk)duk, upon summing over all a. 
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Therefore, using again the representation of /^^{z) as a sum of determinants, Heine's 
formula leads to 

KPn{t-, Z) 

- 1 E n (i - ^) 



= — dct ( fiij - -fJ^ij+i 

T~n \ Z 



0<i,j<n-l 



Z 

= —det{f,ij{t-[z-']))^^..^^_^ 
invoking the fact that 

_ 1 

z 

Formula (3.1.10) follows from computing on the one hand S{t)'mooX{^) using the 
explicit moments /lij, together with (3.1.12), and on the other hand the equivalent 
expression S^-\t)xlz-^). Indeed, using {S{t)xiz))„ ^ p„it; z) = ^oPnkit)z'', 

j>o j>o e>o 



^ ^Z ^^Pniit) / U^+^Pt{u)du 
= j YPniit)u^Y{pj Pt{u)du 



Je 



Ie z — u 
Mimicking computation (3.1.12), one shows 



j>0 
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from which (3.1.10) follows, upon using Snioo = 5^^^. Details of this and subsequent 
derivation can be found in ||. 

The vectors p andq are eigenvectors of L. Indeed, remembering x{.z) = (1, z, z^, ■■■Y , 
and the shift (At>)„ = we have 

J^x{.z) = zx{z) a.Yid Iv' x{.z'^) = zx{.z'^) - zei, with d = (1, 0, 0, ...)^. 

Therefore, p{z) = Sx{z) and q{z) = S^^^xi^^^) ^-^e eigenvectors, in the sense 

Lp = SAS^^Sxiz) = zSxiz) = zp 
L'q = S'^-^X'S'^S'^-\{z-^) = zS'^-\{z-^) - zS^-^ei = zq - zS^-^Ci. 

Then, using L = , one is lead to 

((L — zl)p)n = 0, for n > and ((L — zl)q)n = 0, for n > 1. 

Toda lattice and Lie algebra splitting: The Lie algebra splitting of semi-infinite 
matrices and the corresponding projections (used in (3.1.8)), denoted by ( )g and ( )(, 
are defined as follows: 



gl{oo) =s© b 



3 = {skew-sjTiimetric matrices} 
b = {lower-triangular matrices} 



Conjugating the shift matrix A by S{t) yields a matrix 
L{t) = S{t)AS{t)-^ 

= SAmooS~^ , using (3.1.3), 

= SmooA~^ , using Amoo = mooA^, 

= S{S-^S^-'^)A^S^, using again (3.1.3), 

= {SAs-y = L{ty, 

which is symmetric and thus tridiagonal. Moreover, from (3.1.3) one computes 

. s (a'™. - 

= SA'^S-^ - S^(S-^S^-^)S^ 
otk dtk 

Upon taking the ()_ and ()o parts of this equation {A_ means the lower-triangular part 
of the matrix A, including the diagonal and Aq the diagonal part) leads to 
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Upon observing that for any symmetric matrix 



a c 
c b 



a 
2c b 



a c 
c b 



a c 
c b 



it follows that the matrices L{t), S(t) and the vector p{t] z) = {pnit; z))n>o = S(t)x{z) 
satisfy the (commuting) differential equations and the eigenvalue problem 



f)Q 1 

— = --(L''),S, L{t)p{t-z) = zp{t-z) 



(3.1.13) 



and thus 



dL_ 

dtk 



(Standard Toda lattice) 

The bilinear identity: The functions Tn{t) satisfy the following identity, for n > 
m + 1, t,t' E C, where one integrates along a small circle about oo, 



,it - [^-i])wi(t' + [z-'])e^^'^-'^>'z--"'-'dz = 0. 



(3.1.14) 



An elementary proof can be given by expressing the left hand side of (3.1.14), in 
terms of p„(t; z) and Pm(t, z), using (3.1.9) and (3.1.10). One uses below the following 
identity (see p|) 



f{z)g{z)dz = ( /, 



9[u) 



■du 



z — u 



(3.1.15) 



involving the residue pairing^. So, modulo terms depending on t and t' only, the left 
hand side of (3.1.14) equals 



/ dzz-y^(t-,z)e^^^'^-'>\—'z^+' [ Pui^^Ert'.u^ 

Jz=oo Jr Z U 



p{u)du 



Pn{t] z)e 



Pm{t']z)e^ p{z)dz, using (3.1.15), 



Pn{t; z)pm{t'; z)e^^*^'' p{z)dz = 0, when m < n — 1. 



(3.1.16) 



The residue pairing about z ^ oo between / = J2i>o'^i^'^ ^ ^-'^'^ -9 ^ X^jez^j-^ ^ '= ^ is 
defined as: 



if, 9)00 = f f{z)g{z)— = ^a,h 
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The KP-hierarchy: Setting n = m + 1, shifting t\-^t — y,t'\—>-t + y, evaluating the 
residue and Taylor expanding in i/k and using the Schur polynomials Sn, leads to (see 
footnote 14 for the definition of p{dt)f o g.) 

: 



1 

2Td 


y dz 


1 


j) dz 







)T„ o r„ 



thus yielding (3.1.6), taking into acount the fact that {d / dtijror = and the coefficient 
of yk is trivial for k = 1,2. 

The Riemann-Hilbert problem: Observe that, as a function of z, the integral 
(3.1.10) has a jump accross the real axis 



J_ f Pn(t,M) ^^^^^^^ ^ ^^^^^ z)pt{z) + J— lim 



lz'<0 



2 — *2 

Xz'>0 



Pn{t;u] 



z — u 



pt{u)du, 



and thus we have: (see ^9], |T9|, 0) 
Corollary 3.2 The matrix 



/ r„{t-[z-^) 
' r„(t) ^ 



r„(t) ^ * 



-z 



Tn{t) 



satisfies the Riemann-Hilbert problerr^: 

1. Yn{z) holomorphic on the C+ and 

1 27iipt{z) 

2. Y„(z) = Y„Jz 



(Jump condition) 



3. y:^(z 





19 



C+ and C_ denote the Siegel upper and lower half plane. 



46 



§3, p. 47 

3.2 PfafF Lattice and symmetric/symplectic matrix integrals 

3.2.1 Pfaff lattice, factorization of skew-symmetric matrices and skew- 
orthogonal polynomials 

Consider an inner-product, with a skew-symmetric weight p{y,z), 

{f,9)t= [ [ f{y)9{z)e^'^^y'^''^p{y,z)dydz, with ~p{z,y) = -p{y,z). (3.2.1) 
Since (/, g)t = —{g, f)t, the moment matrix, depending on t = (ti, t2, ■ . .), 

rrinit) = {pij{t))o<ij<n-l = {{y\z^)t)o<i,j<n-l 

is skew- symmetric. It is clear from formula (3.2.1) that the semi-infinite matrix rrioo 
evolves in t according to the commuting vector fields: 



~ f^i+k,j ~r pij-\-kj I.e., 



dt 



dt, 



(3.2.2) 



Since moo is skew-symmetric, rrioo does not admit a Borel factorization in the stan- 
dard sense, but moo admits a unique factorization, with an inserted semi-infinite, skew- 
symmetric matrix J, with = —I, of the form (1.1.12): (see f^) 



moc{t) = Q-\t)JQ'-\t) 



where 
















Q2n,2n 




Q{t) = 




Q2n,2n 








* 


Q2n+2,2n+2 








Q2n+2,2n+2 




I 







G K. 



(3.2.3) 



K is the group of lower-triangular invertible matrices of the form above, with Lie algebra 
t. Consider the Lie algebra splitting, given by 

. , ( t = jlower-triangular matrices of the form (3.2.3)| 
I n = sp{oo) = {a such that Ja J = a|. 
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(3.2.4) 



with unique decomposition^ 



(a_ - J{a+Y J) + -(ao - J^aoY J) 



+ ( (a+ + J{a+yj) + -(ao + J(ao)^J) 



Consider as a special skew- symmetric weight (3.2.1): (see |T3 



P{y, z) = 2D^5{y - z)p{y)p{z) , with a = ^1, p{y) = e'^^^^, 
together with the associated inner- product^] of type (3.2.1): 



(3.2.5) 



(3.2.6) 



{f.9)t = I I f{y)g{z)e^'^^y'^''hD^5{y - z)p{y)p{z)dydz 

f{y)g{z)e^^ *'(^"+'")e(y - z)p{y)p{z)dy dz, for a = -1 

{f,9}{y)e^^''''^p{yfdy, for a = +1 



(3.2.7) 



in terms of the Wronskian {f,g} '■= — /||. The moments with regard to these 
inner-products (with that precise definition of time t\) satisfy the differential equations 
dfiij/dtk = Pi+k,j + yUij+fe, as in (3.2.2). 

It is well known that the determinant of an odd skew-symmetric matrix equals 0, 
whereas the determinant of an even skew-symmetric matrix is the square of a polynomial 
in the entries, the Pfajjianf^ Define now the "Pfaffian r-functions", defined with regard 

'^^a± refers to projection onto strictly upper (strictly lower) triangular matrices, with all 2 x 2 
diagonal blocks equal zero, ao refers to projection onto the "diagonal" , consisting of 2 x 2 blocks. 



^^£(x) = sign X, having the property e' = 2S{x). 

^^with a sign specified below. So det(m2n_i(i)) — and 



(detTO2„(i))'/' = pf{m2n(t)) 



— {dxQ A dxi A ... A dx2n-i) ^ ( fJ-ij{t)dxi A dxj | . 

^' yo<j<j<2n-l J 
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to the inner-products (3.2.7) : 



Pf 



Setting 



p{y)p{z)dydz 



,a = —l 



Pf{ {y\y'}e 



p iy)dy 



0<k,e<2n-l 

a = +l 



0<kl<2n-l 



(3.2.8) 



p{z) = p{z)Ie{z) for a = —1 

p{z) = p^/\z)Ie{z), t t/2 for a = +1, 



in the identities (3.2.8) leads to the identities (3.2.9) between integrals and Pfaffians, 
spelled out in Theorem 3.3 below. Remember c = (ci, ...,C2r) stands for the boundary 
points of the disjoint union _E C M. Sn{E) and %i{E) denotes the set of symmet- 
ric/symplectic ensemble with spectrum in E. 

Theorem 3.3 (Adler-Horozov-van Moerbeke 0, Adler-van Moerbeke 0) The follow- 
ing integrals In(t,c) are Pfajfians: 



« n 

In = / \A^iz)fl\(e^^'^^lpiz,)dz,' 

^Tri-v(x)+zrt.x^)dX, (n even) 



n\pf 

n\Tn(t,C 



£2 



y''z^e{y - z)e^^ ^''■^'^''^ p{y)p{z)dydz 



0<kl<n-l 



(3.2.9) 



/3 = 4 



f eTri-vix)+ETux^)dX, (n arbitrary) 
Jr2„(E) 



nlpf (^j^{y\y^}e^^'^'^^ p{y)dy 

n\T2n{t/2,c). 

The In and Tn 's satisfy 

(i) The Virasoro constraint^ (2.1.7) for = 1,4, 

2r 



0<k,e<2n-l 



- E^^V(q)^ + E (a. ^Jfi,„ - k '^iU^,n) ] Init,c) = 0. (3.2.10) 



^■^ here the a^'s and 6i's are defined in the usual way, in terms of p{z); namely, — 
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(ii) The Pfaff-KP hierarchy: (see footnote I4 for notation) 



Sfe+4(9) - 



1 92 



2 dtidtk 



+3 



r„ O T„ = Sk{d) Tn+2 O T„_2 



(3.2.11) 



n even, /c = 0, 1, 2, .. 



of which the first equation reads (n even) 

Q2 



d 



d 



dtj ^ ^ [dt2) ^dtidts ) ' " V^^? 



logT„ + 6 



92 



12 



Tn-2Tn+2 



(iii) r/ie Pfaff Lattice: The time- dependent matrix L{t), zero above the first superdiag- 



onal, obtained by dressing up A, 



O 



L{t) = Q{t)AQ{t)-' = 



/ * 1 

* {h2/hoy/' 
* 1 

* {hjh2y/^ 
* 



\ 



(3.2.12) 



satisfies the Hamiltonian commuting equations 



dL 



= [-{U)i, L\. (PfafF lattice) 

Uti 



(3.2.13) 



(iv) Skew- orthogonal polynomials: The vector of time- dependent polynomials q{t] z) :- 
{qn{t; z))n>o — Q{'t)x{^) ^ satisfy the eigenvalue problem 



L{t)q{t,z) = zq{t,z) 



(3.2.14) 
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and enjoy the following representations (with = ^'^"^'^It)^ ) 



q2n{t; z) 



''-271 

T2n{t) 



Pf 













Z 














[-1 








2n, -l/2^2n(^M^J) _ ^2n> -1/2 



r2n(t) 



g'2n+i(^; z) 



h 



-1/2 
2n 



r2n(t) 



pf 



( 




1 


/^0,2n-|-l 


m2n{t) 




^2n-l 


/^l,2n-|-l 
/^2n-l,2n+l 


-1 


_^2n-l 





__22n+l 


\ /^2n+l,0 • • • 


Ai2n+l,2n-l 


^2n+l 






Z' 



(3.2.15) 



r/ie?/ are skew- orthogonal polynomials in z; i.e., 

{qi{t;z),qj{t; z))t = Jij. 



The hierarchy (3.2.11) aheady appears in the work of Kac-van de Leur |^ in the 
context of, what they call the DKP-hierarchy, and interesting further work has been 



done by van de Leur |j7TI] 



3.2.2 Sketch of Proof 

Skew-orthogonal polynomials and the Pfaff Lattice: The equalities (3.2.9) be- 
tween the Pfafhans and the matrix integrals are based on two identities the first 
one due to de Bruyn, 

1 f ^ 

\\dy, det(F,(yi) G,{y{) ... F,(|/„) G,(y„)) 
det^/^ ( / iGi{y)F,{y) - F,{y)G ,{y))dy 



0<j,i<2n-l 



and (Mehta Q) 

A^(X) = det (Xl (Xl)' 4 (4)'... (4)')o<.<2n-l 
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On the one hand, (see Mehta [^), setting in the calculation below pt{z) = p{z)e^ Ie 
and 



Fi{x) := / y'pt{y)dy and Gi{x) := F-{x) = x'pt{x), 

J — oo 

one computes: {pt{z) := p{z)e^^^^') 



/ \A2n{z)\Y\pt{Zi)dz., 



i=l 

2n 



/ (^i+l/^*(^i+l))o<i.,<2n-in^^- 

J -00<Zl<Z2<---<Z2n<'Xl ' j^—l 

n 

Y\_Ptiz2k)dZ2k 

00<Z2<Z4,<---<Z2n<00 f^_-^ 



Z2 



det ( / zlpt{Zi)dZi , 4 / 4n-lPt(^2n-l)c?^2n-l , 4n 

'Z2n-2 / 0<i<2n-l 



n 



Y\_Ptiz2k)dZ2k 

00<Z2<Z4<---<Z2n<00 f^_-^ 

det {Fi{z2) , 4 . ^i(^4) - i^i(2;2) , 4 . • • • . i^i(2;2n) - i^i(2;2n-2) , 4n) 



/n 
JJrfzj detfFi(z2) , Gi{z2) , ^^(2:2,,) , Gi{z2r, 

-CO<Z2<Z4<--<Z2n<<Xl ^ 

1 /" " 

- / l[dy, det(F,(yi) , , F,(y„) , 

det^/^ f /" {G,{y)F^{y) - F,{y)G,{y))dy 



0<j<2n-l 



0<i<2n.-l 



0<i,j<2n-l 



= Pfi y'z'e{y-z)e^ruiy^+^np^y)p^,)dydz) = r2„(t), 

establishing the first equation of (3.2.9), taking into account the results in section 2.1. 
On the other hand, upon setting. 



Fj{x) = x^ 



p(x)e^*'"' and G',(x) := F;(x) = (xV(x)e^ , 
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one computes 



1 r " 

l<»J<n k=\ 



0<i<2n-l 



0<i<2n-l 



0<ij<2n-l 



fc=l 

det (xi (xl)' 4 (4)' . . . < (<)') 

1 /" " 

\Je J 0<fe,^<2n-l 

establishing the second equation (3.2.9). 

The skew-orthogonahty of the polynomials ^'^(i; z) follows immediately from the 
skew-Borel decomposition of m^: 

{Qk{t:y),qe{t,z))k,£>o = Q{{y\ z^))i,j>oQ^ = QmooQ^ = J. 

with the QnS admitting the representation (3.2.15) in terms of the moments. 

Using L = QAQ~^, — Q~^JQ^~^ and = — /, one computes from the differ- 
ential equations (3.2.2) 

= (gA^g-i) J - ( jgT-^A^^g^ J) J + ^q-'j - (jq-'^^j)j 

Otk Otk 

T 



Then computing the +, - and the diagonal part (in the sense of (3.2.4) and (3.2.5)) of 
the expression leads to commuting Hamiltonian differential equations for Q, and thus 
for L and q{t;z), confirming (3.2.13): 

^ = -iL%Q, ^ = [{1%, L], ^ = -{L%q. (Pfaff lattice) (3.2.16) 

Oti Oti Oti 

The bilineeir identities: For all n,m > 0, the T2n's satisfy the following bilinear 
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identity 



J z=oo 



+ / + [z])rUt' - [z])e^^<~'^^^~'z'-'"^-^ = 0. 



(3.2.17) 



The differential equation (3.2.2) on the moment matrix rrioo admits the following solu- 
tion, which upon using the Borel decomposition moo = Q^^JQ^^^, leads to: 



(3.2.18) 



and so the right hand side of (3.2.11) is independent of t; say, equal to the same 
expression with t replaced by t'. Upon rearrangement, one finds 



and therefore^ 

(Q(t)x(^) ® {JQ{t')V-'x{z-')) e^Tit.-t',)^^l^ 



27iiz 



{{JQ{t)Y-\{z) ® Q{t')x{z-')) eS?'(*'.-*.)^-'=^. (3.2.19) 



Setting t — t' = [z^ ] + [z2 ] into the exponential leads to 



.):T{t',-tu)--' = ( 1 _ (i- — 



ZZi / \ ZZ2 

and somewhat enlarging the integration circle about 2; = oo to include the points zi and 
Z2 , the integrand on the left hand side has poles at 2; = 2:1 and Z2 , whereas the integrand 
on the right hand side is holomorphic. Combining the identity obtained and the one, 
with Z2 y 00, one finds a functional relation involving a function (p{t; 2;) = 1 + 0{z~^): 

ip{t - [z^^]; zi) ^ ip{t - [z^^]; Z2) ^ ^ (^00 ^ ^ 
(p{t;zi) Vit;z2) 



*using Ax(z) = zx{z), A^x(z) = z ^x{z) and the following matrix identities (see p5[|) 
UiV, - / U,x{z) ® V,^xiz~')T^, U2V2 = i U2x{z) ® v^xiz-') '^^ 



z=oo 27riz J^^o 27riz 
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Such an identity leads, by a standard argument (see e.g. the appendix in [Q) to the 
existence of a function r(t) such that 



T{t) ■ 



This fact combined with the bihnear identity (3.2.19) leads to the bilinear identity 
(3.2.17). 

The Pfaff-KP-hierarchy: Shifting t t-^ t — y,t' ^-^ t + y in (3.2.17), evaluating the 
residue and Taylor expanding in leads to: 



1 
27ri 



+ 7^ / eS?'2,..-V2„+2(t - y + \z\)r2m{t + 2/ - M)^'"-''"rf^ 
2vr2 7^=0 

1 /" °° s °° 



2n-2m-2 



dz 



i=0 

oo 



A:=0 
oo 



— I V^-^s,.(2y)e^-^'^ V^'=s,(9)r2„+2or2, 



5Z Si(-2l/)e^ ^'^'*'Sfe(-9)r2„or2m+2 

j-k=-2n+2m+l 

+ ^ s,(2|/)e^^^'^Sfc(9)r2„+2or2^ 

A:-j=-2n+2m-l 

= ... + yk (^Q^^ ~ Sfe+i((9)^ r2„ o r2„ + Sfc_3((9)r2„+2 o r2„_2^ + 

establishing the PfafF-KP hierarchy (3.2.11), different from the usual KP hierarchy, 
because of the presence of a right hand side. 

Remark: L admits the following representation in terms of r, much in the style of 
(3.1.7), 



h}'\ 





/ ^ 

' -^00 
















Lu 


L12 









* 




L22 


L23 






* 


* 


L32 


L33 
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with the 2x2 entries Lij and /i, being a zero matrix, except for 2 x 2 matrices along 
the diagonah 

For example (' = 



Lr 



-(logrsn)- 1 



_ \ f _ f * (logr2„+2)- 



3.3 2d-Toda lattice and coupled Hermitian matrix integrals 

3.3.1 2d-Toda lattice, factorization of moment matrices and bi-orthogonal 
polynomials 

Consider the inner-product, 

{f,9ks = [ [ f{y)g{z)e^^^'^y'-'^^'^+^y^dydz, (3.3.1) 

on a subset E = Ei x E2 := U^=i[c2i-i, C2i] x Uf^i[c2i-i, C2i] C Fi x F2 C M^. Define the 
customary moment matrix, depending on t = (ti, t2, • • • s = {si, S2, ■ ■ ■), 

1TT'n(t,s) = {fiij{t, s))Q<ij<n-l = ( ) t,s)o<i,j<n-l 

and its factorization in lower- times upper-triangular matrices 

m^it,s) = S^\t,s)S2it,s). (3.3.2) 
Then moo evolves in t, s according to the equations 

Otk OSk Otk OSk 

dM in the integral (3.3.4) denotes properly normalized Haar measure on 7i„. 
Theorem 3.4 (Adler-van Moerbeke 0, |^) The integrals In{t, s; c, c), with Iq = 1, 

r„ = det m„ = i-/„ = -f[ A„(x) A„(y) TT eS?°(*^^'^-^'^^)+^^'=2''=dx,rf|/fc 

•^^^ k=l 

^ ff ^cTriM,M2)^TrETiUMl-s.Mi)^^^^M^^ (334) 
J J HUE) 

satisfy: 
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(i) Virasoro constraint^ (2.2.8) for k > —1, 
1=1 * / 

1=1 ' / 



(3.3.5) 



with 



41 = ^(Jf + (2n + A: + l)J«+n(n + l)Jr), 
J^l = ^(^'^ + (2n + A; + l)jf+n(n + l)Jr). 



(ii) A Wronskian identity^: 

d^logTn d^logTn] , f 9^ log r„ 9^ log r„ 



dtids2 ' dtidsi J , 1 dsidt2 ' dtidsi , ^ 



0. (3.3.6) 



(iii) The 2d- Toda lattice : Given the factorization (3.3.2), the matrices Li := ^lAS*]^ ^ and 
, with hn = ^^r^, have the form (i.e., to be read as follows: the [k — i)th 
subdiagonal is given by the diagonal matrix in front of h^~^) 

= V dtag { !l(^)l!^±^z^±i^^ a^-^ 



00 



h^'h-' = Y.d^ag {^-^t^^^^I^±hzl±l^\ a'^-^, (3.3.7) 



e=o \ rn+k-i+lTn ^^^^ 



and satisfy the 2d-Toda Lattice^ 

^ = [(L?)+,L,] and ^ = [(L^)_ , L,] , z = 1,2, (3.3.8) 



^^For the Hirota symbol, see footnote 14. The J^*'''s are as in remark 1 at the end of Theorem 2.1, 
for /3 = 1 and = 4*'|t— with 



i-\-j—k —i-\-j—k —i~j—k 

^^in terms of the Wronskian {f,g}t — %g — /ff ■ 

^^P+ and P- denote the upper (including diagonal) and strictly lower triangular parts of the matrix 
P, respectively. 
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(iv) Bi- orthogonal polynomials: The expressions 

p^:\t,s;y) := {S^{t, s)x)y)y^ = y 



Tn{t, S) 



(3.3.9) 



form a system of monic bi-orthogonal polynomials in z: 

{Pn\t^ y),Pmit, S; Z))t^s = Sn,mhn With kn = (3.3.10) 

which also are eigenvectors of Li and L2: 

zp^^\t,s;z)=L,{t,s)p^^\t,s;z) and zp^^\t, s; z) = L^it, s)p^^\t, s; z). (3.3.11) 



Remark: Notice that every statement can be dualized, upon using the duahty t < — > —s, 
Li i — > hL^h'^. 

3.3.2 Sketch of proof 

Identity (3.3.4) follows from the fact that the product of the two Vandermonde appear- 
ing in the integral (3.3.4) can be expressed as sum of determinants: 

A„(..)A„(.;) = J2 (3.3.12) 

cr&Sn 



together with the Harish-Chandra, Itzykson and Zuber formula ||3^, ^ 

, ^ n(n — 1) 

/ rf^e-^-^^^" = ^'"^;^^^^^^rr^^ff"- (3.3.13) 
Ju(n) n\ An{x)An{y) 

Moreover the r„'s satisfy the following bilinear identities, for all integer m,n > and 

r„(t - [z-'l .) Wi(t' + [z-'l s')e^^^'^-'>'z''~^''-'dz 



rn+i{t, s - [z\)T^{t\ s' + [z\)e^^^''^-<^'-\^-^~^dz. (3.3.14) 

2=0 

Again, the bi-orthogonal nature (3.3.10) of the polynomials (3.3.9) is tantamount to 
the Borel decomposition, written in the form Simao{hS2~^Y — ^- These polynomials 
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satisfy the eigenvalue problem (3.3.11) and evolve in t,s according to the differential 
equations 

^ = - ((A-'L,A)T") pW ^ = ((fc-iL,fc)T") pW. (3.3.15) 

Otn dSn 

From the representation (3.3.7) and from the bilinear identity (3.3.14), it follows 
that 

Pk-l{dt)Tn+2 O Tn (9^ /o o i «N 
2 = f,^ logTn+i, 3.3.16) 

and so, for /c = 1, 

^--^logr„+i. (3.3.17) 
T^+i dsidti 

Thus, using (3.3.7), (3.3.16) and (3.3.17), we have 



V ^ /n,n+l — — d"^ log r^+i ^ ^ 



Tn+2T, 
Pk-l{ -ds)rn +2 O Tn 

Combining (3.3.18) with (3.3.17) for A; = 2 yields 

•sid 



<9^ logrn+i 
dsidtk 

d"^ logTn + l 

ds\dt\ 
d"^ logrn+i 



dsidti 



(^;L.„« - ^f^^ - 7^>°^f-^) (3.3.19) 



Then, subtracting d/dsi of (3.3.19) from d/dti of the dual of the same equation (see 
remark at the end of Theorem 3.4) leads to (3.3.6). ■ 

3.4 The Toeplitz Lattice and Unitary matrix integrals 

3.4.1 Toeplitz lattice, factorization of moment matrices cind bi-orthogonal 
polynomials 

Consider the inner-product 

{m,g{z))t,s := £ ^J{z)9{z-')e^^^'^^'-^^^-'\ t, s e C^, (3.4.1) 
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where the integral is taken over the unit circle S*^ C C around the origin. It has the 
property 



{z''f,g)t,s = {f,z ^g)t,s. 
The t, s dependent semi- infinite moment matrix moo{t, s), wherep^ 

p{z)dz 



(3.4.2) 



1 27riz 



= Toeplitz matrix 
satisfies the same differential equations, as in (3.3.3): 



0<kl<n-l 



drric 



A^mnr, and 



drrio 
ds„ 



m^k^'\ (2-Toda Lattice) 



(3.4.3) 



(3.4.4) 



As before, define 

Tn{t,s) := detm„(t,s). 

Also, consider the factorization moo(t, s) = S{^(t, s)S2it, s), as in (3.3.2), from which 
one defines Li := SiAS^'^ and L2 := S'2A^5'2'^ and the bi-orthogonal polynomials 
pf'\t,s;z) for k = 1,2. Since moo satisfies the same equations (3.3.3), the matrices 
Li and L2 satisfy the 2-Toda lattice equations; the Toeplitz nature of rrioo implies a 
peculiar "rank 2"-structure, with = 1 — XiVi and Xq = Vq = 1 : 

\ 



and 









1 - xiyi 












-X2yo 


-X2yi 


1 - X2y2 





h-^Lih 






-xsyi 


-X3y2 


1 - x^y^, 










-xm 


-xm 






\ 










( - 




-xoy2 


-xoys 


-Xol/4 




1 - xiui 




-xm 


-xiyA 


L2 = 







1 - X2y2 


-X2y-i 


-X2I/4 












1 - x^y^ 





V 



(3.4.5) 



Some of the ideas in the next theorem are inspired by the work of Hisakado |]37 |. 
A matrix is Toeplitz, when its (i, j)th entry depends on i — j. 
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Theorem 3.5 (Adler-van Moerbeke R]) The integrals In{t,s), with Iq = 1, 



Tn{t, s) = detm„ = : = ^ 
nl n 



U{n) 

J2 sxit)s,i-s), (3.4.6) 

{Young diagrams X \ Ai<n} 



satisfy: 

(i) a SL(2,'L)-algehra of three Virasoro constraints (2.3.2): 

JSC. «) - i%,i-'>. «) - * n) + (1 - 9)J-i,„(-f. ")) s) = 0, 

/or < A; = 0, 6' arbitrary 

y k = i, e = i 

(ii) 2d-Toda identities : The matrices Li and L2, defined above, satisfy the 2-Toda lat- 

tice equations (3.3.8); in particular: 

, Tn-l'Tn+l 

— — -logr„ = 

OSiOti 

and 

\ogTn = -2— log^^ . log r„ - logr„, (3.4.8) 



0S20ti osi r„_i osiOti osfoti 

the first being equivalent to the discrete sinh-Gordon equation 

(In gg„-g„_i _ g<?„+i-<7„^ ^^^^^ ^ log -^"+^ 



(iii) The Toeplitz lattice : The 2-Toda lattice solution is a very special one, namely the 
matrices Li and L2 have a "rank 2" structure, given by (3.4.5), whose Xn's and 
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Un 's equalQ: 

= - [ s„(-TrM,--TrM2,--TrM3,...)e^^^(*>*"'-^»*')t/M 

Tn Juin) 2 3 



U{n) 

a Id Id 



S (-^ )t it S 



Tn{t,s) 

yn{t,s) 

'^n Ju{n) 2 3 

s (JL IJL IJL W it s') 
= ^ =pi^)(M;0), (3.4.9) 

and satisfy the following integrable Hamiltonian system 

dt, "^"^"^ dy^ dt, "^"^"^ dx^ 

^=a-^nyn)^ ^ = -a - ^nVn)^, (3.4.10) 

dSi dyn OSi dXn 

(Toeplitz lattice) 

with initial condition x„(0, 0) = ?/„(0,0) = for n > 1 and boundary condition 
Xo(t, s) = yo(t, s) = 1. The traces 

Ht'' = --TrLl z = 1,2,3,..., fc = l,2. 

t 

of the matrices in (3.4-5) are integrals in involution with regard to the sym- 
plectic structure u := X](r(-'- ~ ^kyk)~^dxk A dyk- The Toeplitz nature ofrrioo leads 
to identities between r 's, the simplest being ( due to Hisakado / p^ ) : 



(3.4.11) 



Remark: The first equation in the hierarchy above reads: 

9Xn I ^ (^Vn \ 

= Xn+l[l - Xnyn) = -|/n-l(l " Xnyn) 

dxji . , Qyn N 

= Xn-\\y - Xnyn) = -yn+l{l - Xnyn) ■ 



^^Remember s(ti,t2, •■•) are elementary Schur polynomials. 
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3.4.2 Sketch of Proof 

The identity (3.4.6) between the determinant and the moment matrix uses again the 
Vandermonde identity (3.3.12), 




Using z^"^ = z ^ (see (3.4.2)), one shows that the polynomials Pn+ii^) ~ zpn\z) 
and pi+i(0)2;"]9n ^2;^^) are perpendicular to the monomials z^,z^, ...,z"' and that they 
have the same 2;°-term; one makes a similar argument, by dualizing 1 <-> 2. Therefore, 
we have the Hisakado identities between the following polynomials: 

pSlW-rf(^) = P^nlMzy^n\z-'). (3.4.12) 

The rank 2 structure (3.4.5) of Li and L2, with Xn — Pn\t, s; 0) and y„ = Pn\t, s; 0), 
is obtained by taking the inner-product of p^^i(z) — zpn\z) with itself, for different n 
and m, and using the fact that zpn\z) — LiPn\z). 

To check the first equation in the hierarchy (see remark at the end of theorem 3.5), 
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consider, from (3.4.9), 

dxn dp'h\t,s;z) 



dti 



dti 



2=0 



2=0 



/lnPn-il(^,S;0)^ 



using (3.3.15), 



i=0 



hi 



using (3.4.5), 



n-l 



7 •^iyi 
i=0 * 



hi hi-i 



using 



hi 

hi-i 



1 •^iUii 



h. 



— Xn+IT 

"-n-l 

•^n+l(l ■^nUn)-! 

and similarly for the other coordinates. From (3.3.7) and (3.4.5), upon making the 
products of the corresponding diagonal entries of Li and hL^h^^, one finds (3.4.11): 



— log ^ log = -Xn+iynXnVn+l 

dti Tn OSi Tn 



h 



- 1- 



hn-1 



1 - 



h 



n+l 



4 Ensembles of finite random matrices 

4.1 PDE's defined by the probabilities in Hermitian, symmet- 
ric and symplectic random ensembles 

As used earlier, the disjoint union E — [c2i-i, C2i] C M, and the weight p{z) — e~^^^\ 
with —p'/p — V — g/f define an algebra of differential operators, {k e Z) 
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The aim of this section is to find PDE's for the following probabilities in terms of the 
boundary points q of E (see (1.1.9), (1.1.11) and (1.1.18)), i.e. 



Pn{E) 



Pn{ all spectral points of M E E) 



n„{E), Sr^iE) or T„{E) 



-tr V(M)^M 



5„( 



n 



k=i 



-tr V(M)^M 



(4.1.1) 



P = 2, 1, 4 respectively. 



involving the classical weights below. In anticipation, the equations obtained in Theo- 
rems 4.1, 4.2 and 4.3 are closely related to three of the six Painleve differential equations: 



weight 




Painleve 


Gauss 




IV 


Laguerre 




V 


Jacobi 




VI 



For (3 = 2, the probabilities satisfy partial differential equations in the boundary points 
of E, whereas in the case /3 = 1, 4, the equations are inductive. Namely, for (3 = 1 (resp. 
/? = 4) , the probabilities Pn+2 (resp. Pn+i) are given in tems of Pn-2 (resp. Pn-i) and 
a differential operator acting on P„. The weights above involve the parameters P, a, b 
and 



S'^ —2 

"1,4: ^ 




-l/2^ 



for p = 2 

1 for /3 = 1,4. 



As a consequence of the duality (2.1.12) between /?-Virasoro generators under the map 
(3 I— >■ 4//3, and the equations (2.1.7), the PDE's obtained have a remarkable property: 
the coefficients Q and Qi of the PDE's are functions in the variables n, (3, a, b, having 
the invariance property under the map 



n 



— 2n, a 



a b 
-2' ^^-2' 



to be precise. 



^ / ^ ^ 0' b. 
Q,{-2n, (3, --,--) 



13=1 



Qi{n,p,a,b)\p__ 



■A ■ 



(4.1.2) 



The results in this section are mainly due to Adler-Shiota-van Moerbeke |TT[ for 
(3 = 2 and Adler-van Moerbeke for (3 = 1,4. For more detailed references, see the 
end of section 4.2. 



65 



§4, p. 66 



4.1.1 Gaussian Hermitian, symmetric and symplectic ensembles 

Given the disjoint union E and the weight e~^^^ , the differential operators Bk take on 
the form 



2r 



d 



Define the invariant polynomials 



Q = m'n (^n + 1 - 1^ and Q2 = 4(1 + dl^)b (^2n + dl^{l - . 



Theorem 4.1 The following probabilities (jS = 2, lA) 

satisfy the PDE's (F := F„ = logPj; 

P„ 2 P„ , 2 \ r 2 when n even and (3=1 



(4.1.3) 



P? 



with index 



1 1 when n arbitrary and (3 = A 
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B\ + {Q2 + m\F)B\ + 4(2 - <4)^(3i3o' - AB.^B^ + 6i3o) ) P. (4.1.4) 



4.1.2 Laguerre Hermitian, symmetric and symplectic ensembles 
Given the disjoint union E C M"*" and the weight z°'e~^^, the Bk take on the form 



2r 



Define the polynomials, also respecting the duality (4.1.2), 

3 
4 

3 

k 2 



Q = 



n(n-l)(n + 2a)(n + 2a + l), for /3 = 1 
n(2n+ l)(2n + a)(2n + a- 1), for /3 = 4 



Q2 



(^(3n^ - ^ + 6an + 4(1 - ^)a + 3^ + (1 - a2)(l - Sl^) 
[(5n'' + 2an + (1 - ^)a) , Qo = h{2 - 5f,J(n + ^). 
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Theorem 4.2 The following probabilities 



satisfy the PDE ^ .■ := F„ = log Pn) 

= {b\ - 2(5f,4 + l)B\ + (Q2 + QB\F - 4(5f 4 + l)B^^F)B\ - ?,5l^{Qi - B^,F)B^, 
+^(2 - 6l,){3Bl - 4SiS_i - 2Bi) + Qoi2BoB., - Bo))f. (4.1.6) 



4.1.3 Jacobi Hermitian, symmetric and symplectic ensembles 

In terms of E G [—1, +1] and the Jacobi weight (1 — the differential operators 
Bk take on the form 

2r „ 



Introduce the following variables, which themselves have the invariance property (4.1.2) 
{Bq = a — b, bi = a + b): 

r = ^{bl + {b, + 2-m s = ^bo{b, + 2-(3) 
4 

qn = -j^iPn + bi + 2- p){pn + 61), 
and the following invariant polynomials in q,r,s: 

Q = — ((s^ -qr + q^f - 4(rs^ - Aqs'^ - 4s^ + gV)) 
16 

= 3^2 _ 3^^ _ ^ + 23g + 24 
Q2 = 3gs^ + - 4g V + 2gr + 4g^ + lOg^ 
Q3 = Sgs^ + 65^ - 3gV + + 4^^ 

Q4 = 9s2_3^^_g^^^2^22g + 24 = gi + (6s2-g2_g). 
Theorem 4.3 T/ie following probabilities 

p .p^_ |An(.)|^ nLl(l - + ^kfdZU 

^ W|A„(^)|/5nLi(i-^;^)"(i + ^^)''rf^^ ^"^ 

satisfy the PDE (F = F^ = logPj: 
^°with the same convention on the indices n ± 2 and n ± 1, as in (4.1.4) 
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for (3 = 2; 



(2Bti + iq-r + A)B\ - {AB^iF - s)B-i + SqB^ - 2qBo + 8BoB\ 

-4(g - l)BiB-i + (4i3_iF - s)Bi + 2{AB-iF - s)BoB.i + 2^-62) F 
+AB\F {2BoF + ^B\F) = 



(4.1.9) 



for (3 = I A: 



Q 




= {q + l) (4(jBii + 12(4^-iF - s)B\ + 2{q+ 12) (4B_iF - s)BqB-i 

Viq^Bl - 4 (g - 4) + q{m-iF - s)Bi + 2{)qBnBl^ + 2q^B^ F 

+ {Q2BI^ - sQiB-i + Q3B0) F + A8{B-iF)^ - 48s{B-iF)^ + 2Q^{B-iFf 
+12 q\BoFf + 16 5 (2 « - 1) {B\F){B^F) + 24 (9 - 1) q{B\Ff 
+24(2B_iF - s) ((g + 2)BoF + (g + 3)B?.iF)b_iF. 



(4.1.10) 



The Proof of these three theorems will be sketched in subsections 4.3, 4.4 and 4.5. 
4.2 ODE's, when E has one boundary point 

Assume the set E consists of one boundary point c = x, besides the boundary of the 
full range; thus, setting respectively E = [—oo,x],E = [0,x],E = [—l,x] in the PDE's 
(4.1.4), (4.1.6) and (4.1.9), (4.1.10), leads to the equations in x below. Notice that, for 
P — 2, the equations obtained are ODE's and, for (3 — 1,4, these equations express 
Pn+2 in terms of P„_2 and a differential operator acting on P„: 

(1) Gauss ensemble {(5 — 2, 1,4) : /„(x) = ^logP„(maxj Aj < x) satisfies 



(2) Laguerre ensemble {(5 — 2,1,4): fn{x) — Xj- logP„(maxj Aj < x) (with all eigen- 




(4.2.1) 



r + 6/;? + (4^(5f,4 - 2) + g^) /; - 4^«, - 2)/, 



n- 
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values Xi > 0) satisfies: 

Sl,Q (^ ^^'l^f' - l] - {^^Jn - ^^{Si, - 2) - Qox - 3Sl,Q^^ U 

- X (^4{6l, + l)fn - - 2) - 2go^ -Q, + 2Sl, + ?j fi 

(4.2.2) 

(3) Jacobi ensemble: f :— fn{x) = (1 — x^)^ logP„(maxj Aj < x) (with all eigenvalues 
— 1 < Aj < 1) satisfies: 

• for P^2 : 

2{x^ - Iff" + 4(a;' - 1) {xf - 3f^) + {l6xf - q{x^ - 1) - 2sx - r) f 

-f{4f-qx-s)^0, (4.2.3) 

• for 13=1,4 : 

= 4{q + l){x^ - 1)2 - 1)/"' + (12/ -qx- 3s)/" + 6q{q - 1)/'^) 

-(a:2 - l)/'(24/(g + 3)(2/ - s) + 8fq{5q - l)x - q{q + l){qx^ + 2sx + 8) + Q2) 
+/(48/3 + 48/2 (ga; + 2x - s) + 2/ {d,q^x^ + 2gx2 - 12gsa; - 24sa; + Q4) 

- q{q + l)x{^qx^ + sx - 2qx - 3q) + Q3X - Qis^ . (4.2.4) 

For P = 2, the term containing the ratio p^" ^ — 1 on the left hand side of (4.2.1), 
(4.2.2) and (4.2.4) vanishes and one thus obtains the ODE's: 

• Gauss: fn{x) := ^ logP„(maxi Aj < x) satisfies: 

+ 6 /'2 + 46(2n - bx^)f' + 4b^x / = 

• Laguerre: fn{x) := log P„(maxj Aj < x) satisfies 

x^f" + xf + 6xf - Aff - ((a - hxf - Anhx)f - b{2n + a - bx)f = 0. 
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Jacobi: = (1 — x^)^ log P„(maxi Aj < x) satisfies: 



2{x^ - Iff" + A{x^ - 1) {xf - 3/'2) + (I6x/ - qix^ - 1) - 2sx - r) f 

-f (4/ -qx-s) = 0. 

Each of these three equations is of the Chazy form (see section 9) 

r + ^f" + |r - %ff + + - ^/ + ^ = 0, (4.2.5) 

with c — and P,Q,R having the form: 

Gauss P{x) = 1 4Q(x) = -462^2 + 86n i? = 

Laguerre P{x) = x 4Q(x) = — (6x — a)^ + 46na; i? = 

Jacobi P{x)^l-x^ 4:Q{x) ^ -^{q{x'^ - 1) + 2sx + r) R^O 

Cosgrove shows such a third order equation (4.2.5) in f{x), with P{x), Q{x), R{x) 
of respective degrees 3, 2, 1, has a first integral (9.0.2), which is second order in / and 
quadratic in /", with an integration constant c. Equation (9.0.2) is a master Painleve 
equation, containing the 6 Painleve equations. If f\x) satisfies the equations above, 
then the new (renormalized) function g{z), defined below. 

Gauss g{z) = b-^'"^ f{zh-^l'^) + Inz 
Laguerre g{z) = f{z) + |(2n + a)z + ^ 
Jacobi g{z) := -y{x)\^=2z-i - f + 

satisfies the canonical equations, which then can be transformed into the standard 
Painleve equations; these canonical equations are respectively: 

• /2 ^ _4^/3 ^ ^^g> -gf^ A^g' + A2 (Painleve IV) 

. {zg"f = {zg' - g) [-Ag'^ + A,{zg' - g) + A^) + + A, 

(Painleve V) 

. {z{z - l)g" f = (zg' - g) [Ag'^ - Ag\zg' - g) + A^) + A,g'^ + A,g' + A, 
with respective coefficients (Painleve VI) 

. A^ = b\ A2 = b\{n + 1)2 + 4), A3 = -a%{n + f ), = ^ 

•((^ + §)' + *) 
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For {3 = 1 and 4, the inductive partial differential equations (4.1.4), (4.1.6), (4.1.10), 
and the derived differential equations (4.2.1), (4.2.2) and (4.2.4) are due to Adler-van 
Moerbeke 0. For (3 = 2 and for general E, they were first computed by Adler-Shiota- 
van Moerbeke []TT|, using the method of the present paper. For j3 = 2 and for E 
having one boundary point, the equations obtained here coincide with the ones first 
obtained by Tracy- Widom in ^4 , who recognized them to be Painleve IV and V for the 



Gaussian and Laguerre distribution respectively. In his Louvain doctoral dissertation, 
J. P. Semengue, together with L. Haine |3^, were lead to Painleve VI for the Jacobi 



ensemble, for (3 = 2 and E having one boundary point, upon subtracting the Tracy- 
Widom differential equation (@]) from the one computed with the Adler-Shiota-van 
Moerbeke method (|TT|]). Cosgrove's ([^) and Cosgrove-Scoufis's classification (||24|, 



(A. 3),) leads directly to these results. 

4.3 Proof of Theorems 4.1, 4.2 and 4.3 
4.3.1 Gaussian and Laguerre ensembles 

The three first Virasoro equations, as in (2.1.29) and (2.1.32), are differential equations, 
involving partials in t G C°° and partials Vi,V2,'D3 in c = (ci,...,C2r) G M^*", for 
E := En(t, c) = login', they have the general form: 

dF 

'^^^=7^+ Yl IkjVjiE) + + Sktu A; = 1,2, 3, (4.3.1) 

-i<j<k 



with first Vj(F)'s given by: 

'dii+j ' 2""'" \ dt'i ' \dtij 



i,i+j>l 



In (4.3.1) and (4.3.2), (3 > 0,7^^,7^,5^ are arbitrary parameters; also 62j = for j 7^ 2 
and = 1 for j = 2. The claim is that the equations (4.3.1) enable one to express all 
partial derivatives. 



d'^+-+ikF{t, c) 



dt\\..dtl' 



, along C := {all ti = 0, c= (ci, ...,C2r) arbitrary}, (4.3.3) 

c 



uniquely in terms of polynomials in 

Indeed, the method consists of expressing dE / dtk\^^Q in terms of Pfc/I^^g, using (4.3.1). 
Second derivatives are obtained by acting on T>kE with T>i, by noting that commutes 
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with all ^-derivatives, by using the equation for P^F, and by setting in the end t — 0: 

dF 



VeVkF = Ve 



dtk 



-l<j<k 



m I ^H-P"). provided Vj{F) does not 



-i<j<fe 



contain non-linear terms 



dF 



7)7:+ E Yl l^^V,iF) + Seh 



dtk 

Vf 

dtkdte 



-l<j<k 



-i<j<e 



+ lower- weight terms. 



When the non-linear term is present, it is taken care as follows: 



dF OF 
2—Ve— 
dti dh 



dF d OF 



-i<j<e 



Higher derivatives are obtained in the same way. We only record here, for future use, 
the few partials appearing in the KP equation (3.1.6): 



d^F 



d^F 



dtj 
d^F 



dtl 



= (X>^ - 7ioI?i) + 7io7i - 

= {V\ - 67102?? + ll7i'o^? - 6710^1) F - 67i'o(5i - 71710) 

V% - 272oI'2 -I- /372i732l'f - ((271 -I- 710)721732/? -I- 272,-i)I'i 



d'^F 



dtidts 



-2721P3 j F + P^2i732{'DiFf 

-|-/?72i732(7i + 7io7i - ^1) + 2(72173 + 72072 + 7i72,-i) 



T^iT^a - ^732©? -I- /?732(7i + 2710)©? - ^710732(271 -|- 710)^1 



-37i,_ir>2 - 37ioI?3 ]F+ Y7io732(2?ii^)^ - I3'^s2{T>iF){VIF) 
+^(271073 + /?7327io(7i + 7io7i - + 271,-172). 



4.3.2 Jacobi ensemble 

Here, from the Virasoro constraints (2.1.35), one proceeds in the same way as before, 
by forming BiF\^^Q , BiBjF\^^Q, etc., in terms of ti partials. For example, from the 
expressions , , BqF\^^q , one extracts 



OF 




d^F 




OF 








dtj 


t=o 


^2 


t=o 
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Prom the expressions B^^FL , BqB. 




=0 ' '-^1-' lt=0' 


tion, one extracts 








dF 






d^F 




t=o 


dtl 





t=0 
2 ; 



Finally, from the expressions B2F\^^f^ , BiB-iF\^^q , BqF\^^q , BoBiiF\^^^ , BliF\^^^, 
one deduces 



d^F 




OF 








d^F 




d^F 


dtf 


t=o 




t=o 


dtldt2 


t=o 


dtidts 


t=o 


dtl 



(4.3.4) 



t=o 



This provides all the partials, appearing in the KP equation (3.1.6). 
4.3.3 Inserting pcirtials into the integrable equation 

Prom Theorem 3.3, the integrals In{t,c) , depending on P = 2,1,4, on t = {ti,t2,---) 
and on the boundary points c — (ci, C2r) of E, relate to r-functions, as follows: 



In{t, C) 



|A„(^) 



k=l 



n\Tn{t,c), n arbitrary, 13 — 2 
n\Tn{t, c), n even, P — ^ 
n\T2n{t/2, c), n arbitrary, (5 — 4 



(4.3.5) 



where Tn(t,c) satisfies the KP-like equation 



,o T„-2(t,c)r„+2(^,c) ^ J n arbitrary for /3 = 2 

12 ()i^4 = log r„(t,c). 



with 



T„(i, C)2 



d 



d 



n even for j3 — 1,4 



(4.3.6) 



dtl J \dt2 / UL\UL^ 

Evaluating the left hand side of (4-3.6): Here Init) will refer to the integral (4.3.5) over 
the full range. Por [3 = 2, the left hand side is zero. Por (3 = 1, the left hand side 
can be evaluated in terms of the probability Pn{E), as follows: taking into account 

Pn:=Pn{E)=In{Q,c)/In{Q), 



12 



rn-2(t,c)Tn+2{t,c) 



= 12- 



'ri-2 



(t,c)In+2(t,c) 



t=0 



t=0 



12 



n-2)!(n + 2)! 4(t,c)2 

n(n-l) /„_2(0)/„+2(0)P„_2Pn+2 



(n + l)(n + 2) 4(0)2 
,.u(i) Pn-2{E)P^+2{E) 



p2 

n 
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with bn^ given byQ 



n{n-l) /„_2(0)/n+2(0) 
(n + 2)(n+l) 4(0)2 



n(n— 1) 
1662 



(Gaussian) 



n(n-l)(?i+2a)(n+2a+l) 
16P 

(Laguerre) 
(Jacobi) 



(4.3.7) 



For /5 = 4, we have: 



12 



T-2n-2(t/2,c)r2„+2(i/2,c) 



r2„(t/2,c)^ 



t=0 



12 
12 



[nl 



/„_i(t,c)/„+i(t, c) 



(n-l)!(r2 + l)! 4(t,c)2 
n /„_i(0)/„+i(0)P„_iP„+i 



(n+1) 4(0)2 p2 
(4)iVi(^)iW^ 

p2(^) ' 



with 



/n-l(0)4+l(0) 



1)1(^ + 1)! /2(o) 



^"4^2^^^ (Gauss) 

2ra(2n+l)(2n+a)(2n+a~l) 

(Laguerre) 
(Jacobi) 



(4.3.8) 

where Q is precisely the expression appearing on the left hand side of (4.1.10), and 



where Qq is given by 



Qf = 3q(g + l)(9-3)(g + 4±4v/^ 



+ for P = 1 
- for /3 = 4 



(4.3.9) 



The exact formulae bn^ and bn' show they satisfy the duality property (4.1.2): 
bl^\a,b, n) = 6«(-|, ~ -2n). 



^"'^this calculation is based on Selberg's integrals (see Mehta p 340). For instance, in the Jacobi 
case, one uses 



W r(a + J/3/2 + l)r(6 + J/3/2 + l)r((j + 1)13/2 + 1) 
r(/3/2 + l)r(a + 6+(n + j-l)/3/2 + 2) 
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Evaluating the right hand side of (4-3.6): From section 2.4, it also follows that Fn(t; c) = 
\ogIn(t; c) satisfies Virasoro constraints, corresponding precisely to the situation (4.3.1), 
with 

Gaussian ensemhl^^ 

71 - 1 = 7i,o = 7i = 0, (5i = -f 

72- 1 = 0, 72,0 = -1/2, 72,1 = 0, 72 = -f 0-1, ^2 = 

73,_1 = -^CTi, 73,0 = 0, 73,1 = -i, 73,2 = 73 = 0, (^3 = -f CTl. 

Laguerre ensembl^. 6i = 62 = 6^ = 0, and 

7i,_i = 0, 7i,o = -1, 7i = -f (0-1 + a), 
72-1 = 0,72,0 = -CTl, 72,1 = -1,72 = -f 0-1 (0-1 + a), 
73,_i = 0,73,0 = -0-10-2,73,1 = -0-2,73,2 = -1,73 = -f 0-10-2(0-1 + a). 

Jacobi ensemble: see (4.3.4). 
They lead to expressions for 











d^F 




d^F 


dtf 


t=0 


dtl 


i=0 


dtidt^ 


i=0 


dtl 



in terms of T>k and Bk-, which substituted in the right hand side of (4.3.6) - i.e. in the 
KP-expressions - leads to the right hand side of (4. 1.4), (4. 1.6), (4.1.9) and (4.1.10). In 
the Jacobi case, the right hand side of (4.3.6) contains the same coefficient l/Qg as in 
(4.3.9), which therefore cancels with the one appearing on the left hand side; see the 
expression h]i'^ in (4.3.7) and (4.3.8). 



5 Ensembles of infinite random matrices: Fredholm 
determinants, as r-functions of the KdV equation 

Infinite Hermitian matrix ensembles typically relate to the Korteweg-de Vries hierarchy, 
itself a reduction of the KP hierarchy; a brief sketch will be necessary. The KP-hierarchy 
is given by t^-deformations of a pseudo-differential operator^ L: (commuting vector 
fields) 

dL 8 

— = [(L")+,L], L = D + a_,D-^ + with D = —. (5.0.1) 



■^^ Remember from section 2.1, ai = (3{n — 1 
^■^Remember from section 2.1, ai = 
^''in this section, given P a pseudo-differ 
(strictly) smoothing part of P respectively. 



^■^Remember from section 2.1, cri = — 1) + a + 2 and (T2 = f3{n — |) + a + 3. 
^''in this section, given P a pseudo-differential operator, P+ and P_ denote the differential and the 
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Wave and adjoint wave functions are eigenf unctions t; z) and \l/~(a;, t; z), depend- 

ing on X G M, t G C°°, 2; G C, behaving asymptotically like (5.0.3) below and satisfying: 

zm+ = LVI/+, -- = (L")+M/+, 

J|l (5.0.2) 
;^vl/- = L^vl/-, __ = -(lT")+x[/-. 

According to Sato's theory, \1'+ and \E'~ have the following representation in terms of a 
r-function (see [p5| ): 



T[t) 

= e^(""+^?°*'^")(l + 0(^~')), for 2/00, (5.0.3) 
where r satisfies and is characterized by the following bilinear relation 

/,Er.-...V,<-|.-.lM..|.-,)..^0. fo.a..,.eC~; (5.0.4) 

the integral is taken over a small circle around z = oo. From the bilinear relation, one 
derives the KP-hierarchy, already mentioned in Theorem 3.1, of which the first equation 
reads as in (3.1.6). 

We consider the p-reduced KP hierarchy, i.e., the reduction to pseudo- differential 
Us such that U = + . . . is a differential operator for some fixed p > 2. Then 
(L^^)+ = for all A; > 1 and thus dL/dt^p = 0, in view of the deformation equations 
(5.0.1) on L. Therefore the variables tp,t2p,t3p, ... are not active and can thus be set 
= 0. The case p = 2 is particularly interesting and leads to the KdV equation, upon 
setting all even tj = 0. 

For the time being, take the integer p > 2 arbitrary. The arbitrary linear combina- 
tions0 

$^(x, t; z) := a^^^(x, t; uz) (5.0.5) 

are the most general solution of the spectral problems = and L^^^^ = 

respectively, leading to the definition of the kernels: 

/X 
dx ^~{x,t;y)^^{x,t;z), and k^^{y,z) := k^^t{y,z)lE{z), (5.0.6) 

where the integral is taken from a fixed, but arbitrary origin in R. In the same way 
that \E'^(x,t, z) has a r-function representation, so also does k^^{y,z) have a similar 
representation, involving the vertex operator: 

Y{x,t;y,z) := ^ a-a^,X{x,t;ujy,uj'z), (5.0.7) 



35 



:= {w such that = 1} 



76 



^5, p. 77 



1 a 



(5.0.8) 



where (see ||2^, [TT|, |T! 

A condition ^^^g^^ ^^^77^ = is needed to guarantee that the right hand side of (5.0.7) is 
free of singularities in the positive quadrant {yi > and Zj > with i, j = 1, n} and 
Ymiy^z t; y, z) exists. Indeed, using Fay identities and higher degree Fay identities, 
one shows stepwise the following three statements, the last one being a statement about 
a Fredholm determinant^: 



r{ty 



l<j,j<n 
E 



1 

-Y\Y{x,t;yi,Zi)T 

T 

i=l 



det(/ - Aej 



~A dzY {x,t;z,z) 



T 



T 



r{t,E) 
T{t) ■ 



(5.0.9) 



The kernel (5.0.12) at t = will define the ststistics of a random Hermitian ensem- 
ble, when the size n 00. The next theorem is precisely a statement about Fredholm 
determinants of kernels of the form (5.0.12); it will be identified at t = with the prob- 
ability that no eigenvalue belongs to a subset E; see section 1.2. The initial condition 
that Virasoro annihilates tq, as in sections 2.1.2 (Proof of Theorem 2.1), is now replaced 
by the initial condition (5.0.11) below. 

r(2) 



Theorem 5.1 (Adler-Shiota-van Moerbeke |TT|, [T^) Consider Virasoro generators J\ 
satisfying 



d_ 

dz 



1 .(2) 



jr{t),Y{x,t-z,z) 



(5.0.10) 



where Y{x,t; z, z) is defined in (5.0.7), and a r-function satisfying the Virasoro con- 
straint, with an arbitrary constant Ckp: 



(^■Jkp ~ Cfcp) r = for a fixed k > -1. 
Then, given the disjoint union E C M^, the Fredholm determinant of 



<t(A,A'):=-%^/i.(A'), X = z^,X' 
P z 2 z 2 



(5.0.11) 



(5.0.12) 



^^The Fredholm determinant of a kernel A{y, z) is defined by 

det(/ -\A) = l+f2 (-^)'" /• • ■/ det (^(^^, ^j))i<,,,<„, dzi . . . dz„ 



2i <■ ■ ■<2:„ 
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satisfies the following constraint for that same k > —1: 



(- E + ^(4? - j r det(J - f^K^,) = 0. (5.0.13) 



The generators JA take on the following precise form: 

4+j = Tl 



E a5^ + 2 E ^*^|-+ E (5.0.14) 



Remark: For KdV (i.e., p = 2), {L'^)~^ = = — g(x) holds, and thus the adjoint 
wave function has the simple expression: \E'~(a;, t; z) = '^~^{x,t; —z). In the next two 
examples, which deal with KdV, set 

^'(x,t; z) := ^+(x,t; z). 

Example 1: Eigenvalues of large random Hermitian matrices near the "soft 
edge" and the Airy kernel 

Remember from section 1.2, the spectrum of the Gaussian Hermitian matrix ensemble 
has, for large size n, its edge at ±^/2n, near which the scaling is given by ^J2n}l^ . 
Therefore, the eigenvalues in Theorem 5.2 must be expressed in that new scaling. Define 
the disjoint union E = U^[c2j_i, C2i], with C2r possibly oo. 

Theorem 5.2 Given the spectrum zi > Z2 > ... of the large random Hermitian matrix 
M, define the "eigenvalues" in the new scale: 

= 2ni(^^-l) forn/oo. (5.0.15) 

The probability of the "eigenvalues" 

P{E^) := P{ all eigenvalues Ui G E^) (5.0.16) 
satisfies the partial differential equation {setting Bk := Yl^i=i(^^^^^)Tl 

[b\ - 4(^0 - \)^ S-i \ogP{E') + Q{B\ \ogP{E-)f = 0. (5.0.17) 
^''When C2r = cxD, that term in Bk is absent. 
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In particular, the statistics of the largest "eigenvalue" ui (in the new scale) is given by 

P{ui < x) = exp j {a-x)g^{a)da^ , (5.0.18) 

with 

g" = xg + 2g^ (Painleve II) 

-2j . (5.0.19) 

^^(^) - -t/h^ forx / oo. 

The partial differential equation (5.0.17) is due to Adler-Shiota-van Moerbeke |[Tl|, 
T2[. The equation (5.0.19) for the largest eigenvalue is a special case of (5.0.17), but 



was first derived by Tracy- Widom [Q, by methods of functional analysis. 

Proof: Remember from section 1.2, the statistics of the eigenvalues is governed by the 
Fredholm determinant of the kernel (1.2.4), for the Hermite polynomials. In the limit, 

1 / / — u t — V \ 

lim — — -K^ V2n + — — -, V2n + — — - = K(u, v), 

where 

/•oo /"OO 

K{u,v)= A{x + u)A{x + v)dx, A{u) = I e^"^-^'/^cia;. (5.0.20) 

Jo J -oo 

Then 

P(E'=) := P( all eigenvalue Ui E E"") = det(/ - K{u, v)Ie{v)). (5.0.21) 

In order to compute the PDE's of this expression, with regard to the endpoints q of 
the disjoint union one proceeds as follows: 

Consider the KdV wave function \E'(x, t; z), as in (5.0.2), with initial condition: 

v|>(x, to; z) = z\A[x + z^") = e^^+f ^'(1 + Oiz-')) , ^ ^ oo, to = (0, 0, ^, 0, ...), 

(5.0.22) 

in terms of the Airy function^, which, by stationary phase, has the asymptotics: 

-| poo 3 3 

A(m):=— / e-V+y«rfi/ = M-2ei"^(l + 0(M-t)). 



The definition of A{u) is slightly changed, compared to (5.0.20). A{u) satisfies the 
differential equation A{y)" = yA{y), and thus the wave function \E'(x,to;-2) satisfies 
{D^ - x)^{x,to;z) = z^^{x,to;z). Therefore L'^\t=to = SD^S-^\t=tQ = - x, so 



38 



The i in the definition of the Airy function is omitted here. 
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that is a differential operator, and \1/ is a KdV wave function, with r(t) satisfyin,_|__| 
the Virasoro constraints (5.0.11) with C2k = The argument to prove these 

constraints is based on the fact that the hnear span (a point in an infinite-dimensional 
Grassmannian) 



f ^2 3 d"- 

W = spanc <^ ^niz) := e 3^ ^"^'^^^ 



n = 0,1,2, 



is invariant under multiphcation by z"^ and the operator + 2z'^) — jz ^. 
Define for A = A' = ^'^ the kernel KtiX, A'), 

1 r°° 

Kt{X,X') := , ^ '${x,t;z)'^{x,t;z')dx, (5.0.23) 
2^:2 2;' 2 Jo 

which flows off the Airy kernel, by (5.0.22), 



1 



Kt, (A, A') = 2 / + + 



'0 



Thus r det(/ — i^T^J satisfles (5.0.13), with that same constant Ckp, for A; = —1, 0, 1, ... : 
E^r^ A + 1 + ^det(/ - K^) = 0. (5.0.24) 



Upon shifting ts 1— > + 2/3, in view of (5.0.22), the two flrst Virasoro constraints 
for A; = -1 and A; = read: {8^ := E?Ii c- +'^) 

^-.logr(,i.) = (|^ + ^E^t.^)log.(t,i^)4 

The same method as in section 4 enables one to express all the t-partials, appearing 
in the KdV equation, 

■^^ Although not used here, the r-function is Kontsevich's integral: ((see |T^) 

T(t) = ^ with tn = Tr{Z ") + -(5„ 3, Z = diagonal matrix. 
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in terms of c-partials, which upon substitution leads to the partial differential equation 
{B^^ - 4(i3o - i)) / + 6(i3_i/)2 = (announced in (5.0.17)) for 

/ := logPiE'^) = E ^ log^(^'), where P(i?=) = det(/ - K^^) = 
When E = (— oo, x), this PDE reduces to an ODE: 

- Axf + 2/ + 6/'^ = 0, with / = A log P(max A, < x). (5.0.26) 

ax i 

According to Appendix on Chazy classes (section 9) , this equation can be reduced to 

+ 4/'(/'' - xf + /) = 0, (Painleve II) (5.0.27) 

which can be solved by setting 

/' = and f = g"- xg' - g\ 

An easy computation shows g satisfies the equation g" = 2g^ + xg (Painleve II), thus 
leading to (5.0.19). 

Example 2: Eigenvalues of large random Laguerre Hermitian matrices near 
the "hard edge" and the Bessel kernel 

Consider the ensemble of n x n random matrices for the Laguerre probability distribu- 
tion, thus corresponding to (1.1.9) with p(dz) = z^^'^e'^^'^dz. Remember from section 
1.2, the density of eigenvalues near the "hard edge" 2 = is given by 4n for very large 
n. At this edge, the kernel (1.2.4) with Laguerre polynomials pn tends to the Bessel 



kernel [52 



lim —K^^;^ (—^—\= K^^\u,v) :=- [ xJJxu)JJxv)dx. (5.0.28) 
n/tx, An \An An/ 2 Jq 

Therefore, the eigenvalues in the theorem below will be expressed in that new scaling. 
Define, as before, the disjoint union E = U^[c2i_i, C2i]. 

Theorem 5.3 Given the spectrum < zi < Z2 < ... of the large random Laguerre- 
distrihuted Hermitian matrix M, define the "eigenvalues" in the new scale: 

Ui = Anzi for n oo. (5.0.29) 

The statistics of the "eigenvalues" 

P{E^) := P( all "eigenvalues" Ui e E") (5.0.30) 
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leads to the following PDE for F = log P{E^): {setting Bk := Eili cf^^^) 

Bt - 2Bl + (1 - u'')Bl + B^[b^-^^F~ A{BoF){BlF) + Q{BlFf = 0. (5.0.31) 

In particular, for very large n, the statistics of the smallest eigenvalue is governed by 

— / du I , ui ~ Anzi, 

Jo u J 

with f satisfying 

{xf'f -A[xf' - f)f'^ + ((x - - /)/ = 0. (Painleve V) (5.0.32) 
The equation (5.0.32) for the smallest eigenvalue, first derived by Tracy- Widom []64| , 



by methods of functional analysis, is a special case of the partial differential equation 



(5.0.31), due to In], p 



Remark: This same theorem would hold for the Jacobi ensemble, near the "hard edges" 
z = ±1. 

Proof: Define a wave function \l/(a;,t; z), flowing off 

^(x, 0; z) = e'^'Bi-xz) = e^"(l + 0{z-^), 
where B{z) is the Bessel function^ 

B{z) = e^ze^HSz) = ^p^;^ / (^I^^Tj^^^ = ^ + Oi^^')- 
As the operator 

L |t=o — 5 

is a differential operator, we are in the KdV situation; again one may assume t2 = t^ = 
• • • = and we have 

T{t) 

in terms of a r-function^ satisfying the following Virasoro constraints 

j£V = ((2z/)2 - l)4or. (5.0.33) 



e = ^yV2e*''''/^ -1/2 <u < 1/2. 



T{t) is given by the Adler-Morozov-Shiota-van Moerbeke double Laplace matrix transform, with 
tn given in a similar way as in footnote 32 (see |lo[] ): 

T{t)^c{t)f dXdetX^-i/2e-Tr(z^^) /■ dr^Q(y)e-Tr(xy^) ^ 
Jnl Jh+ 
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Set p = 2, aj" = all = i^^e'^''^^ and aZi ^ af ^ ^e-'"""/^ in (5.0.5); this defines the 
kernel (5.0.6) and so (5.0.12), which in terms oi X — z'^ and A' = z''^, takes on the form: 

Ki:h\X') = r (^e'^^*{x,t,z) + e-'^^*{x,t,-z))■ 

A^^Vzz' J ^ f 

e~~'^{x,t,z') +ie~'^{x,t, -z') \ dx, 

which flows off the Bessel kernel, 

Ki"!>{X,X') ^ l f xMxV\)J,{xy/^)dx. 
^ Jo 

2(A-A0 ^~ 
The Predholm determinant satisfies for E C IR+ and for A; = 0, 1, . . . : 

- E + + - 5m) {r det(/ - <f ) = 0. (5.0.34) 

Upon shifting i— > + -\/^ and using the same Bi as in (5.0.25), the equations 
for A; = and k — 1 read 



Bologr{t,E) = Uj2^t~ + ^^]^ogr{t,E) + 

\i>i * V 

B.logr(t,£) = -(^|:.(._ + -^ + V^_ + -_j.ogx(t,£). 

(5.0.35) 

Expressing the i-partials (5.0.20), appearing in the KdV-equation at i = (see 
formula below (5.0.25)) in terms of the c-partials applied to logT(0,£'), leads to the 
following PDE for F = logP(£;'^): 

Bt - 2Bl + (1 - u^)Bl + i3i (^i3o - 0^ F - A{BoF){BIF) + Q{BlFf = 0. (5.0.36) 

Specializing this equation to the interval E = (0, x) leads to an ODE for / := —x dF j dx, 
namely 

r + -/" - -r + -,ff + ^^^/' - ^/ = 0, (5.0.37) 

X X x'^ x^ 2x^ 

which is an equation of the type (9.0.1); changing x r\ —x and / rv — / leads again to 
an equation of type (9.0.1), with P{x) — x, 4(5(x) — —x — v'^ and R — Q. According to 
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Cosgrove-Scoufis ||2^ (see the Appendix on Chazy classes), this equation can be reduced 
to the equation (9.0.2), with the same P,Q,R and with c = 0. Since P{x) = x, this 
equation is aheady in one of the canonical forms (9.0.3), which upon changing back x 
and /, leads to 



rr + 4{-xf + f)f" + {{x - v^)f -f)f = 0. (Painleve V) 



[X 

Example 3: Eigenvalues of large random Gaussian Hermitian matrices in 
the bulk and the sine kernel 

Setting u = ±1/2 yields kernels related to the sine kernel: 

^-,0 [y ) - ^ ]^ ^1/2^1/2 

1 /sinx(?/ — z) sm-xiy + z) 



271 \ y — z y + z 

;^(-l/2). 2 ^2N _ 1 r COSXy COSXZ 
^-,0 12/ ' ^ J - ^J^ ^1/2^1/2 

1 f sin x{y ~ z) _^ smx{y + z] 



2tt \ y — z y + z 

Therefore the sine-kernel obtained in the context of the bulk-scaling limit (see 1.2.7) 
is the sum K^^q^'^'^ +kI. q^'^^. Expressing the Fredholm determinant of this sum in terms 
of the Fredholm determinants of each of the parts, leads to the Painleve V equation 
(1.2.8). 



6 Coupled random Hermitian ensembles 

Consider a product ensemble (Mi, M2) G := 7i„ x 7i„ of n x n Hermitean matrices, 
equipped with a Gaussian probability measure, 

CndM,dM2 e-fTr(Mf+M|-2cMiM2)^ g.l) 

where dMidM2 is Haar measure on the product 7i^, with each dMj, 

n n 

dMi = Al{x) Y[dxidU and dMa = Al{y)Y[dyidU (6.0.2) 

1 1 

decomposed into radial and angular parts. In terms of the coupling constant c, appear- 
ing in (6.0.1), and the boundary of the set 

E = EixE2:= U[=i[a2i-i,a2i] x UtiN-i,&2i] C M^ (6.0.3) 
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define differential operators Ak, Bk of "weiglit" k, 

, \^ d d „ -r^, d d 



i=i ^ i=i ^ 



forming a closed Lie algebraQ. The following theorem follows, via similar methods, 
from the Virasoro constraints (3.3.5) and the 2-Toda equation (3.3.6): 

Theorem 6.1 (Gaussian probability) (Adler-van Moerbeke 0) The joint statistics 



Pn{M e nl{E, X E2)) 



nl{EixE2) 



I [ A„(x)A„(y)ne-i(-i+^'-2-'=?^'=)f/xfcd?/, 
k=i 



[[ A„(x)A„(y)TTe- 
J Jr2. ^J, 



satisfies the non-linear third-order partial differential equatior^ (independent of n) 
{Fn:=i\ogPn{E) ): 

B^AiE^ , BiAiFn + ^^\ - (^2SiF„ , A^ii^n + ^^l =0. (6.0.4) 

7 Random permutations 

The purpose of this section is to show that the generating function of the probability 

P{L{T:n) <e) = ^#{7r„ G Sn \ L(7r„) < £} 
nl ' 



42 

1 — 1 — 

[A2,B2] = [Ai,B2] = —^Bi [Bi,B2] = l^Bi. 

*'^in terms of the Wronskian {/, g^x — X f.g — f.Xg, with regard to a first order differential operator 
X. 
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is closely related to a special solution of the Painleve V equation, with peculiar initial 
condition. Remember from section 1.4, L(7r„) is the length of the longest increasing 
sequence in the permutation 7r„. 

Theorem 7.1 For every i > 0, 



n=0 



u{i) Jo ^uJ 



(7.0.1) 



with gi satisfying the initial value problem for Painleve V.- 

9 - -7r[ 7 + - + - + -g{g - 1 - TTi- 

2 \g-l gj u u 2u^ g 



(7.0.2) 



iu) = l 



u 



+ 0(m^+^), nearu = 0. 



The Painleve V equation for this integral was first found by Tracy- Widom ||6^ . This 
systematic derivation and the initial condition are due to Adler-van Moerbeke 

Proof: Upon inserting {ti,t2,.--) and (si,S2,...) variables in the f7(n)-integral (7.0.1), 
the integral 



In{t, s) 



(7.0.3) 



n\ det 



dz 
27riz 



n!r„(t, s) 



0<k,e<n-l 



puts us in the conditions of Theorem 3.5. It deals with semi-infinite matrices Li and 
hLjh~^ of "rank 2", having diagonal elements 



Or,. := TT— lor 



{Li)n-i,n-i and 6* : = 



\0g^^ = {hLlh ^)n-l,n-l- 



Otl Tn-1 dsi r„_i 

To summarize Theorem 3.5, /„(t, s) satisfies the following three types of identities 
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(i) Virasoro: (F := logr„) (see (3.4.7)) 

. V_1T„ 



= 



Vor„ 



i>i 



i>2 



dsi 



VlTn 



i>l 



j>2 



= 



dti Tn 



q2 q2 q2 



i>l ^ " i>2 



dtidsi dtidsi 



(ii) two-Toda: (see (3.4.8)) 



(7.0.4) 



logTn 
ds2dti 



o ^ 1 



g2 ^3 

log - TTJ^ log '^ri 



= 26: 



9si Tn-i dsidti dsfdti 

Q2 Q3 



(7.0.5) 



(iii) Toeplitz: (see (3.4.11)) 

^{Tju = — log — lo. 



dti Tn-1 dsi Tn-1 



+ 1 + 



dsidti 

= -bnbl+(l + 



52 d f d T„ 

logT„ ) ( 1 + - — -— logr„ - — — log 



dsidti 



dsi \dti Tn-1 



(7.0.6) 



Defining the locus C = { aA\ ti — Si — 0, except 7^ 0}, and using the second 

relation (7.0.4), we have on C, 



d a 

c \ ati osi 



0, 



implying r„(t, s) 



is a function of x := — tiSi only. Therefore we may write r„ 



Tn{x), and so, along C, we have 



a 



19^' didll - Settmg 



(7.0.7) 
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and using x — —tiSi, the two-Toda relation (7.0.5) takes on the form 
9Mogr„ 



si- 



ds2dti 



dsidti 



dsi \ dsidti 



.hi 



Setting relation (7.0.7) into the Virasoro relations (7.0.4) yields 



= 







d 



d 



= Uit^-Sit;— log 



dti Tn 



dti dsi 
92 



Tn-l 



= tibn + Sibl 



■Si 



n- 



ds2dti dtidsi 



logr„ 



+ n 



= -x(2^/4x) + /;(x))+n(-/„(x) + l). 



This is a system of two linear relations in 6„ and 6* , whose solution, together with its 
derivatives, are given by: 



n 
tl 



bn_ n(/„-i) + x/; db,, _ d K _ xju: - z;^) + (/n + n)/; 

Si 2x/„ ' dsi dx si 2/2 



Substituting the result into the Toeplitz relation (7.0.6), namely 

d 



bJl - (1 - /„) g^Kj , 

leads to /„ satisfying Painleve equation (7.0.2), with g = fn, as in (7.0.7) and u — x. 
Note, along the locus C, we may set ti — ^/x and si = since it respects si = —x. 

Thus, /„(t, s)|^ equals (7.0.1). 

The initial condition (7.0.2) follows from the fact that as long as < n < £, the 
inequality L(7r„) < £ is always verified, and so 



(n!) 

thus proving Proposition 7.1. 

Remark: Setting 



X 



n\ (£ + 1)!' 



,((£ + l)!-l) + 0(a;^+2) 



exp I X — 



X 



(^+1)!2 



fn{x) 



g{x) 



g{x) - 1 

leads to standard Painleve V, with a — S — 0, (3 — 7 = —2. 
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8 Random involutions 

This section deals with a generating function for the distribution of the length of the 
longest increasing sequence of a fixed-point free random involution 712]^, with the uniform 
distribution: 

P <i + l, < G 5°,) = e si I L«) <i + l}. 

Proposition 8.1 (Adler-van Moerbeke 0) The generating function 



00 



fc=0 

771 ^xTrM I 771 ^xTrM 

= -c^o(^+i)_e +-c/o(£+i)+e 



exp ( /" i^du] + exp f ^l^du 1 , (8.0.1) 



'0 M y \Jq u 

where f = , satisfies the initial value problem: 

r . i/" + i/''^ - ±ff - 15^/ . 15; , ^Al^ . 

=u'^± + 0(m^+2), near m = 0. (Painleve V) 

(8.0.2) 

Proof: The first equality in (8.0.1) follows immediately from proposition 1.1. The 
results of section 1.3 lead to 

I e'^^'^dM = 6^'' [ A„{zff\e^'''^{l-ZkT{l + Zk)''dzk, (8.0.3) 

io(2n+l)± -'[-1,1]" k=l 

with a = ±1/2, b = ^1/2, (with corresponding signs). Inserting tj's in the integral, the 
perturbed integral, with e^^ removed and with ti = 2x, reads 

/n 
A„(^)2 17(1 - ZkTil + Zkfe^^'^'''^dzk = n!r„(t); (8.0.4) 

this is precisely integral (3.1.4) of section 3.1.1 and thus it satisfies the Virasoro con- 
straints (3.1.5), but without boundary contribution BiF. Explicit Virasoro expressions 
appear in (2.1.35), upon setting [3 = 2. Also, T„(t), as in (3.1.4), (see Theorem 3.1) 
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satisfies the KP equation (3.1.6). Differentiating the Virasoro constraints in ti and t2, 
and restricting to the locus 

C :— {ti — X, all other ti = 0}, 
lead to a linear system of five equations, with bo = a — b, bi = a + b, 



1 






(j 


d 


1 






d 


1 




in 



/ In 



■^1+2 



in five unknowns (F„ = logr„) 



+ &oJi+i + J 



(1) A J 

k+2 J 



0, A; = -1,0 



0, A; = -1,0 



0, k^-l 







OF,, 




d^F^ 




d'F^ 




d'F,, 


dt2 


c 




c 


dtidt2 


c 


dtidts 


c 


dtl 



Setting ti — X and F^ — dFn/dx, the solution is given by the following expressions. 



dt2 
dFn 



dtidt2 
d'F„ 



dhdts 



-^({2n + hi)F^ + n{bo-x)) 

--^{x {f;1 + F'^ + {bo - + n{n + bi)) - (2n + 6i) ((2n + 6i)F; + 6on)) 

-^((2n + 6i)(x<'-f;)-6on) 

(x^K" + 2F;f;') - X {{x' - box + + + 6oF; + (2n + b.^F^l 



+n(n + + 2(2n + bifF^ + 2bon{2n + bi) 



^ (2F;:2 + 26o< + ((2n + bif + 2)^;;' + 2n(n + 

-3(2n + 6i)2f; - 36on(2n + . 



Putting these expressions into KP and setting ti — x, one finds: 



d 



- 7T- 



9 



^2 



dtidts 



Fn + G(T^F,, 



9t? 



= ^ ( x^F"" + 4a;2F'" + x {-4x^ + 4boX + 2-{2n + bif) F" + Sx^F'F" 
+6x=^F"^ + 2xF''^ + (26ox - (2n + F' + n{2x - bo){n + 6i) - 60^^)- 
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Finally, the function H{x) := x-^F{x) = x-^\ogTn{x) satisfies 

x'^H'" + xH" + QxH'"^ - [AH + Ax^ - Ahx + {2n + af) H' + (4a; - 2h)H 

+2n{n + a)x — bn{2n + a) = 0. 

(8.0.5) 

This 3rd order equation is Cosgrove's |2^, equation, with P = x, AQ = —Ax"^ + 
Abx — {2n + a)^, 2R = 2n{n + a)x — bn{2n + a). So, this 3rd order equation can 
be transformed into the Painleve V equation (9.0.3) in the appendix. The boundary 
condition /(O) = follows from the definition of H above, whereas, after an elementary, 
but tedious computation, /'(O) = /"(O) = follows from the differential equation (8.0.5) 
and the Aomoto extension (see Mehta [0, p. 340) of Selberg's integrahQ 

. . I' \A{x)f UU ^J(l - x,ydx,...dx„ ~l{j + S + 2+i2n-j- l)/3/2- 

However, the initial condition (8.0.2) is a much stronger statement, again stemming 
from the fact that as long as < n < £, the inequality L(7r„) < £ is trivially verified, 
thus leading to 

(2 i+1 \ 

T^(£TT)!+^^"'"'V' 

ending the proof of Proposition 8.1. ■ 



9 Appendix: Chazy classes 

Most of the differential equations encountered in this survey belong to the general Chazy 
class 

= F{z, /, /', /"), where F is rational in /, /', /" and locally analytic in z, 

subjected to the requirement that the general solution be free of movable branch points; 
the latter is a branch point whose location depends on the integration constants. In his 
classification Chazy found thirteen cases, the first of which is given by 



, 1 i?e 7 + 1 i?e 5 + 1 

''"where Rej, Re5> -1, Re (3 > -2min -, , 

n n — 1 n — I 
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with arbitrary polynomials P{z), Q{z), R{z) of degree 3, 2, 1 respectively. Cosgrove and 
Scoufis ||2^, (A. 3), show that this third order equation has a first integral, which is 
second order in / and quadratic in /", 

+ l{P"f + Q")f - \p"'f + c) = 0; (9.0.2) 
c is the integration constant. Equations of the general form 

r = Gix, /, /') 

are invariant under the map 

aiz + 02 , „ a^f + UqZ + a-j 
X I— > and / (-^ . 

a^z + 04 032; + 04 

Using this map, the polynomial P{z) can be normalized to 

P{z) = z{z — 1), z, or 1. 

In this way, Cosgrove shows (9.0.2) is a master Painleve equation, containing the 6 
Painleve equations. In each of the cases, the canonical equations are respectively: 

• = -Af> - 2g'{zg' - g) + Ai (Painleve II) 

• g"2 = -4^'3 + 4(zg' - gf + A^g' + A2 (Painleve IV) 

. {zg" f = izg' - g) (-^^ + A.izg' - g) + A^) + A,g' + A, 

(Painleve V) 

. iz{z - l)g'r = {zg' - g) {Ag'^ - Ag'{zg' - g) + A2) + A,g'^ + A^g' + A, 

(Painleve VI) 

(9.0.3) 

Painleve II equation above can be solved by setting 

= ^(«+(«^+|)^^l)^ (^ = ±1). 
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Then u{z) satisfies yet another version of the Painleve II equation 

u" = 2u^ + zu + a. (Painleve II) 

Now, each of these Painleve II, IV, V, VI equations can be transformed into the standard 
Painleve equations, which are all differential equations of the form 

/" = ^(z, /,/'), rational in /, /', analytic in z, 

whose general solution has no movable critical points. Painleve showed that this re- 
quirement leads to 50 types of equations, six of which cannot be reduced to known 
equations. 
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